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TORSIONAL VIBRATION 
IN DIESEL ENGINES 


INTRODUCTION 


ToRSIONAL vibration in diesel engines is due primarily 
to the fact that the periodically varying cylinder torque 
acts upon an elastic shaft system which itself has one 
or more inherent natural frequencies, depending upon 
the proportions of the shaft and the number and 
distribution of the attached masses. Each frequency 
is distinguished by the number of nodes formed in the 
shaft during vibration, the lowest frequency being that 
associated with one node, the lowest but one with two 
nodes and so on. 

When the shaft is rotating at such a speed that the 
power strokes of the engine synchronise with one of 
the natural frequencies of the shaft system, it is said 
to be in a critical speed. When the speed is such that 
the natural frequency is exactly twice the number of 
revolutions per second, the shaft is said to be in a 
critical speed of the second order. Thus, in a four cycle 
engine, it is possible to have critical speeds of every 
integral order and half order and in a two cycle engine 


critical speeds of every integral order, the order being 
1 ] 
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determined by the number of vibrations occurring per 
revolution of the engine. 

Critical speeds are not all of equal importance. In 
general, the most important are those the torque 
components of which act cumulatively, so as to twist 
the shaft in its normal mode of vibration. These are 
known as “ major” criticals, all others being termed 
‘“minor”’ criticals. The major critical speeds most 
frequently encountered are those the order of which 
is a direct multiple of the number of power strokes per 
revolution. In the majority of cases, minor critical 
speeds are not important, though in some instances 
they may set up stresses sufficiently high to fracture 
a crankshaft and for this reason they must always be 
taken into consideration. The actual stress set up by 
a critical speed is dependent on a number of factors 
which will be considered later. 


THE TORSIONAL PENDULUM 


Free Harmonic Vibrations.—The torsional pendulum, 
which is the simplest example of a system capable of 
torsional vibration, consists essentially of a disc 
supported by a shaft as shown in Fig. 1. 

If the disc is given a smal] initial twist and then 
suddenly released, torsional vibrations will ensue. 
These vibrations, which are maintained by the elastic 
restoring forces in the shaft, are known as the “ free” 
or “‘ natural ” vibrations of the system and an expression 
for their characteristics can be found by setting down 
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the differential equation of motion of the system and 
obtaining its solution. 

Let [ =the polar moment of inertia of the disc. 
For a uniform disc of radius of gyration 7 ft. and weight 
W lbs. 


2 
I= a Ibs. ft.2 sec.? 


Let R = the stiffness of the shaft. For a uniform 
shaft of length 7 ft. and diameter d ft., the material 


I 


Fig. 1.—Torsional Pendulum. 


of which has a torsional elastic modulus G@ Ibs. per 


sq. ft. 


Grd‘ 
R= — ft. Ibs. per radian. 


Let « = the angle of twist at any instant. Then the 
torque moment in the shaft at any instant will be Re. 
Neglecting the mass of the shaft,* this torque moment 
will be equal in magnitude but opposite in sign to the 

* In nearly all cases, it is possible to neglect the mass of the shaft 
when dealing with torsional vibration in diesel engines, as the effect on 
the frequency is very small. If it is desired to take the shaft into account, 


one third of the moment of inertia of the shaft should be added to that 
of the mass. 
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moment of inertia of the disc times its angular accelera- 
tion, whence 
d?a 


[aa + ha = 0 


or Ie+Ra=0 . ‘ . (1) 


Dividing through by J and letting 2 = w? 


atwa=od0 . : ~ (2) 
The solution of this equation is 


a = Cie + Cre-* where i == V—1 
= A cos wi + B sin wt j i ~- (3) 


where A and B, which equal (C, + C,) and a(C, — C,) 
respectively, are arbitrary constants. 

This solution shows that the vibration of the dise is 
Simple Harmonic Motion. 

By writing 4 =a sin g and B =a cos 9 we have 


a=a sin (wt + ¢) ‘ . (4) 


This is the general equation for Simple Harmonic 
Motion, in which @ is the phase velocity, usually in 
radians per second, a is the amplitude of the oscillation 
and is the phase angle when ¢ = 0. 


: Rk : es ee 
Now since w* = +, the motion of the dise is given by 


pean (#1 + ») 4 () 


and since w is the phase velocity, the time of one 
complete vibration is 

2 
=~ seconds 
a 


t= 


THE TORSIONAL PENDULUM 5 
and the frequency of the vibrations 
1 @ 


t= eo On 


1 ae er 
=sV7 vibrations per second . (6) 





Free Vibrations with Damping. 
considered, the effect of damping has been neglected, 
and it was assumed on page 2 that once the torsional 
pendulum has been disturbed, the vibrations will 
continue indefinitely with constant amplitude. In 
actual practice damping, in the form of hysteresis in 
the material, air resistance, etc., is always present, 
and the natural vibrations gradually decrease in 
amplitude until the pendulum finally comes to rest. 
The actual relationship between this damping force 
and the other factors involved is at present an open 
point, but in dealing with the torsional pendulum, it 
is usually assumed for mathematical convenience that 
the damping is proportional to the velocity of the 
vibrations. 

Therefore if u is the coefficient of friction in foot-lbs. 


. a 
per radian per second, the damping force is Ba 


The differential equation for damped free vibration 


now becomes 
7 Oo 


‘ie +g eee 
or Ie+uce + Ra =0 ; . (7) 
Dividing through by J and letting 


we have 
a++2atwe=O0 . ; . (8) 
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The general solution of this equation is 
ase (C,e'V ev sk Coe! at pit) 
provided that p? is less than w?. If p? is not less than 


w? the motion is dead beat and non-vibratory. 

Therefore 

a = e-(A cos Vw? — pt + Bsin Vw? — p%) (9) 

where A and B, which equal (C, -+- C,) and 1(C, — C,) 
respectively, are arbitrary constants. 

This represents simple harmonic motion whose 
amplitude diminishes asymptotically to zero. 

Equation (9) can be represented as simple sinusoidal 
motion. 

Let A =a sin g and B =a cos @ where a and 
are constants as in equation (4). | 


Then «=e7~"{asin (Vw? — pt + —)} . . (10) 
As in the case of Equation (5) the natural frequency is 
given by 

L, .. poset 
f=sVot—pt.  . . (11) 


If p is small compared with w, p? can be neglected, 
and the frequency then becomes the same as that for 
undamped simple harmonic motion. It may thus be 
assumed that the natural frequency is unaffected by 
the presence of a small damping coefficient. 


Substituting - for w* and S for p we have 


1 /R B\*? , 
f= NT (5 vibrations per second (12) 


The ‘rate of diminution of amplitude” of the 
damped torsional pendulum can be very conveniently 
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expressed in terms of the “logarithmic decrement,” 
Referring to Fig. 2, which shows the amplitude plotted 
against wi, the time between A and B is 


Q7r 
= — seconds 
w 


In this time, therefore, the amplitude drops from 


e~" to er(t+7) 





Fig. 2.—Logarithmic Decrement. 


27p 
The second amplitude is the first multiplied by e * 


and since this factor is the same for any two consecutive 

maxima, the amplitudes decrease in geometric series. 
But numbers in geometric series have their logarithms 

in arithmetic series. Thus the difference between the 


: . 2 
logarithms of two consecutive maxima, 1.¢. a , Is 


constant. This is known as the “ logarithmic 
decrement.” 


' Substituting . for 2p we have 


Logarithmic Decrement = = . . (18) 
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Forced Vibrations with Damping.—Consider now the 
case when the damped simple harmonic motion of 
Equation (7) is acted upon by a periodic disturbing 
force of the form a, sin gi. 

The differential equation of motion now becomes 


d*a da. 
er + he + Ra =a, sin gt 
or Ia + wa + Ra = a, sin gt. . (14) 
Dividing through by I and letting 
Ce a ee ae 
[=~ 2D 7 =e and 7 = 5 


we have 
a+2pa+we=b sin gt . . (15) 


The complcte solution of this equation is 
a =e asin (Vw? — pt + y)} + Psin(qt--e) . (16) 
b 


where P= V(o* — g2)? + 4pg? (17) 
“pq 
and tan é= (a — 9°) ‘ . - (18) 


The first term on the right-hand side of Equation (16) 
represents the natural vibrations, previously considered, 
which are soon damped out. The second term 
represents the forced vibrations. 

When the natural vibrations have been damped out 

a = P sin (qt — «) 

This is simple harmonic motion of the same frequency 
as the disturbing force, and lagging an angle ¢ behind 
it. The maximum amplitude is P. 


Now it will be remembered that b = + and w? =+ 


b a; 
so that ot R 
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a, . : ; 
Now R is the maximum value of the disturbing 


force divided by the shaft stiffness. This ratio thus 
represents the amplitude which would be produced by 
the steady application of the force ay. 
Therefore the maximum “ static ”’ amplitude is 
b 
6 = oa : ‘ . (19) 
Dividing the top and bottom of Equation (17) by 
w®, the maximum “ dynamic ” — becomes 


Gey a V4 7 Va-BY EE (20) 


where f = oi = the natural frequency. 
7 == qg/2x —= the forcing frequency. 
= 2p/w == a constant depending on the damping 
forces. 
The fraction 


(0 ay aE . (21) 


therefore represents the degrec of ade of the 
static amplitude due to the relationship between the 
three factors f, f,; and K. It is commonly known as 
the “dynamic magnifier.” 

If f, is very small compared with f, M = 1 and the 
** dynamic ” amplitude equals the “ static ”’ amplitude, 
t.e. the amplitude may be calculated on the basis that 
the disturbing force is steadily applied. 

If f, =f, M=1/K. Thus if the damping is small, 
the dynamic amplitude is very large and under these 
conditions a small disturbing force may produce an 
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amplitude of sufficient magnitude to fracture the shaft. 
This is the important condition of resonance. 
If f, is very large compared with f, M tends to equal 
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zero, and the dynamic amplitude thus actually becomes 
much smaller than the static amplitude. 

These three conditions are most clearly shown in 
the form of a diagram. Fig. 8 shows the dynamic 
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magnifier M (Equation 21) plotted against the ratio 
filf for various values of K. What may be termed 
the ‘“‘ danger zone ” is clearly secn to lie approximately 
between the ratios f,/f = 0-8 and 1-20. 

As the ratio f,/f varies, the angle of lag e also varies. 
Referring to Equation (18), if q is small compared with 
w, the angle of lag is small. If g = @, the condition for 


780 
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Fig. 4.—Relationship between Angle of Lag ¢ and Ratio f,/f. 


resonance, € = 7/2, t.e. the phase of the vibrations is 
90° behind the disturbing force. If q is large compared 
with w, tan ¢ is negative, t.e. the phase is opposite to 
that of the disturbing force. 

_ Fig. 4 shows these conditions diagrammatically and 
it will be seen that for low values of the damping forces 
there is a very rapid change in phase as resonance is 
approached. 
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THE NATURAL FREQUENCIES OF MULTI-MASS 
SYSTEMS 


The Frequency Equation.—The torsional pendulum 
considered on page 2 forms a basis for all diesel engine 
natural frequency calculations, but only in extreme 
cases such as a single cylinder engine with a heavy 
fly-wheel or a radial aero engine with a heavy airscrew 
is it directly applicable to the solution of torsional 
vibration problems. In all other cases, more than one 
mass and, in most cases, more than one connecting 








Fig. 5.—System Comprising (7 -+ 1) Masses and n Shafts. 


shaft, are involved, so that the extension of the principles 
already considered will now be applied to multi-mass 
systems. 

Referring to Equation (1) it will be seen that the 
natural frequency of the torsional pendulum was found 
by setting down the differential equation of motion of 
the system and obtaining its solution. From this 
equation the natural frequency was found directly 
(Equation 6). 

To turn now to the system shown in Fig. 5. Damping 
will be neglected as experience shows that the amount 
of damping encountered in diesel engines ‘is of such 
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magnitude that the natural frequency is unaffected 
thereby. The system comprises (rn + 1) masses and 
m connecting shafts. J,, Z,, etc., are the moments of 
inertia of the masses and #,, R,, etc., are the shaft 
stiffnesses respectively. Let «a,, a, etc., be the 
instantaneous angles of twist of the masses. 

Then the following equations of motion hold good :— 


T,a, + Ry(a, — a2) = 0 
0, + Bla, — 43) — R(x, — a.) = 0 
I3%3 + Rg(ag — a4) — Bo(a. — as) = 0 


I, at a R, ie ios O41) Soe F(%,-1— a, ) = 0 


non 


Fiyr%nat — R,,(« a, — O41) = 0 (22) 


From this set of differential equations, a general 
equation for the natural frequency can be obtained.* 
By writing Equations (22) in a slightly different form, 
the same equation can be obtained by an alternative 
method.f It is not proposed to repeat these proofs 
here, but merely to state the equation, which is as 
follows :— 


F" — F"-'9, +j2 + Js tja t+ + Jan) 
+ F"~ 25 (j3 +a + e+ Jon) FJalJa + ++ -Jan) toot 
— Fj, t9s(Js - - oo + JalJo-++Jan) + .-++} 
+ Jaljalje + + Jan) 1 oeeee} Hjsl-.-e} Ht ecee] 
SO os oie SS eee aan =O. : - (28) 


where F' represents the natural frequency (squared) of 


* Lochner, J/.J.H.H#., 1927, Vol. 65, p. 76. 
t Bradbury, Proc.J.Mech.H., 1933, Vol. 125, p. 483. 
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the whole system, and j,, j,, etc., are the individual 
frequencies (squared), 1.¢ 


jz is the frequency (squared) of mass m, on shaft 1 


J2 99 99 39 99 Me 99 1 

ds 99 99 99 9 Meo 99 2 

Is 99 99 9 99 Ms 29 2 
etc. 


The equation is built up as follows. The first term 
contains F' to the nth power, where n is the number 
of connecting shafts. The second term contains F to 
the power (rn —1) multiplied by the sum of the 
individual frequencies (squared). The third term con- 
tains F' to the power (n — 2) multiplied by the sum of 
the individual frequencies (squared) up to j2,~ >, each 
of which is first multiplied by the sum from 7,4, to 
Jens aNd SO on. 

The equation has n roots each of which represents 
the square of a natural frequency of the complete 
system. Thus the important fact is established that 
a system such as that shown in Fig. 5 has possible 
natural frequencies equal in number to the connecting 
shafts. 

Now consider the case where n = 1. This gives the 
simple two mass-one shaft system shown in Fig. 6. 
Putting n = 1 in Equation (23) we have 


oa (91 + je) = 0 
whence Fe= ie +J.  . . (24) 


Substituting 7, = (=); j= (2) 


f, 
{AV Ay 
- &) +7] 


NATURAL FREQUENCIES OF MULTI-MASS SYSTEMS 15 


whence the natural frequency 


fedyR(Et+z) ©. on 


The form or “ mode” of vibration of this system is 
shown by the dotted line in Fig. 6. When the mass J, 
is swinging in one direction, the mass J, is swinging 
in the opposite direction, and there is therefore a point 
in the shaft at A where the vibrations change phase 
and are of zero magnitude. This point is known as 
a “ node.” 





Fig. 6.—T wo Mass-One Shaft System. 


The position of the node is found as follows :—The 
system of Fig. 6 is virtually two torsional pendulums 
of equal frequency joined back to back, the point of 
juncture being the node. Letting the total length of 
the shaft be / and the lengths on either side of the 
node J, and I, respectively, we have 

1 /Ri 1 /Rl 
f= 3,V Gi, ~ VTL 
whence Li-=i4; 


I, ly 
ane oe 


16 TORSIONAL VIBRATION IN DIESEL ENGINES 


Thus the shaft is divided in mverse proportion to 
the moments of inertia of the masses. 


EXAMPLE.—-Two Mass-One Shaft System. 
The system which is shown in Fig. 6 has the following 
specification :— 
I, = 86:3 lbs. ft.? sec.? 


I, = 100 lbs. ft.? sec.? 
R= 4-7 x 108 ft. Ibs. per radian. 


t == 10 inches. 
By Equation (24) 
F = Jy “+ Je 
1 \? 4:7 x 108 
n= \an) —ge3 ~ = 3280 
1 \? 4:7 x 108 
a= (=) 499° 


F = 9, + Jo = 4470 
Natural frequency 
f = V4470 = 66-8 v.p.s. 
Or by Equation (25) 
pnt 4-7 x 10° x 186-3 
2m 100 x 36°3 


= 66°8 V.p.s. 
To find the position of the node 
L 4, 
i, = I, = 0 aoe 
whence l, = 0-3681, 
and since 1=1, +1, = 10 ins. 
10 


i= 1-363 7°85 ins, 
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Now consider the three mass-two shaft system shown 
in Fig. 7. Putting n = 2 in Equation (23) we have 


F? — F(j, +Jj2t+Js + ja) + [9s(Js + 5a) + I2Ja] = 0 (26) 


giving two values for F and hence two natural 
frequencies. 

The lower of these two natural frequencies will occur 
with one of the end masses swinging against the other 
two masses, the node lying in one of the two connecting 
shafts. The higher natural frequency will occur with 





Fig. 7.—Three Mass-Two Shaft System. 


the two end masses swinging in phase, the centre mass 
swinging opposite in phase. In this case the nodes 
will lie one in each connecting shaft. 

These two modes of vibration are shown by the 
dotted lines in Fig. 7. 

The modes of vibration again depend on the relative 
magnitudes of the masses and shaft stiffnesses. In 
general it may be assumed that if (j, +.) is less than 
(73 +34), the node for the lower frequency will lie in 
the first shaft and vice versa. Also if I, =I, and 

2 
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R, = R,, the node will lie in the centre mass. The 
determination of the actual positions of the nodes is 
best carried out by the method described in the 
following section on calculation by successive approxima- 
tion. 


ExaMPLE.—Three Mass-Two Shaft System. 
The system which is shown in Fig. 7 has the following 
specification :-— 
I, = 86 lbs. ft.? sec.? 
I, = 36 lbs. ft.? sec.? 
I, = 100 lbs. ft.? sec.? 
R, == 5-02 x 108 ft. Ibs. per radian. 
R, = 5-38 x 108 ft. Ibs. per radian. 


By Equation (26) 
| ale F(j, + Je + Js + Ja) a [ji(Js + ja) + JaJa] = 0 


The insertion of the appropriate values in this equation 
is best carried out in tabular form, as follows :— 


ee ee ee oo nee ee rm ee ED ee 


R 1\? R 
: ~ T (=) T \Ostie)| 44 
dt 36 5:02 « 108 13:9 x 104 3520 
Js 36 5-02 x 108 13-9 x 104 3520 
ds 36 5-38 x 108 14:9 x 104 3770 3770 
N 100 5-38 x 106 5:38 x 104 1360 1360 | 1360 


ATT | ES | AE | Ss | I. | EPS | uthrARNS ETE 


Acasa sae er Soe ten RFOEASMY Sems smart 
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f2 — F(12170) + 3520(5180) + 8520(1860) = 0 
F? — F(12-17 x 108) + 22-8 x 10¢ =0 
which factorises as 
(F — 2320)(F — 9850) = 0 


The frequencies thus are 


One node. 
f; = V2820 = 48-2 v.p.s. 


Two node. _ 
fo = V9850 = 99-3 v.p.s. 


The two modes of vibration, which give the positions 
of the nodes, are shown by the dotted lines. As 
previously mentioned, the determination of the positions 
of these nodes 1s best carried out by the method 
described in the following section. 

Frequencies by Successive Approximation.—Whilst it 
is possible to evaluate the natural frequencies of a 
system comprising any number of masses by means of 
the equation already given, the method is apt to become 
very laborious when n is greater than 2. For this 
reason it is usual, when calculating systems having 
four or more masses to use the method of “ successive 
approximation.” 

Referring again to Equations (22), these equations 
will have solutions of the form 


“, = a, sin (wt + eé,) 
, = A, sin (wt + &), 
etc. 


where w is the phase velocity of the particular frequency 
under consideration. 
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Substituting these values of «,, a, etc., in Equations 
(22) we have: 


L,a,0° => R,(a, =— Ay) —_ 0 
T,a,w* — R,(a, — a3) + R,(a, — a2) — 0 
I,ayo* — R(ags — ay) + R,(a, — as) — 0 


[4,07 — k,(4, — 4,41) + Ry (a, 1— @,) = 0 


an 


I, $1441 a? = R,(a, —~ G4 1) =- 0 (27) 


For calculation purposes, these equations are more 
conveniently written in the following form :— 


“wt 
dg == a, — R, (Z,4;) 


at 
ag = a, — R, (1,4, + 1,a2) 


a 


a, = as — R, (L,a, + La, + I as) 


2 
G..1,=% — = (Iya, + Ia, +...1,a,)  ~ (28) 


Now, by making an estimate of the value of w?, and 
assuming an arbitrary value, say unity, for a,, the 
angular displacement of the first mass, the value of a, 
can be calculated from the first equation. a, having 
been found, a, can be calculated from the second 
equation and so on up to the last equation, which gives 
the value of a,,,, the angular displacement of the 
last mass. If w? has been correctly chosen, the equation 


[,a,0? + I a0? +... + Dn 4 14y 4 10" == Q (29) 
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representing the sum of Equations (27), will be satisfied. 
If this equation is not satisfied, fresh values of w? 
must be tried until satisfaction is obtained. 

The following example shows a simple method of 
setting out this calculation in tabular form, the system 
being that of a six cylinder diesel engine as shown 
in Fig. 8. The specification is as follows :— 


Moment of inertia per cylinder 
line. : é . == 7-11 Ibs. ft.” sec.? 


Moment of inertia of flywheel = 400 Ibs. ft.? sec.* 
Stiffness between cylinder centres = 14-1 X 10° ft. lbs. 


per radian. 
Stiffness between No. 6 cylinder 


and flywheel ; ; . = 14-25 x 10° ft. Ibs. 
per radian. 





Fig. 8.—Six Cylinder Diesel Engine. 


Column 1 gives the moments of inertia of the masses 
in ft. lb. second units, while Column 2 gives the moment 
of inertia multiplied by the square of the angular 
velocity @. Column 8 gives the angles of rotation of 
the masses, that of mass No. 1 being given an arbitrary 
value of unity. Column 4 gives the moments of the 


inertia forces of the masses and Column 5 the total 
(Table overpage) 
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twisting moment of the inertia forces of all the masses 
to the left of the scction under consideration. The 
shaft stiffnesses are given in Column 6 and these, divided 
into the twisting moments given in Column 5, give the 
angles of twist for consecutive parts of the shaft 
(Column 7). 

To satisfy Equation (29) the final figure in Column 5 
should be zero. To obtain this exactitude, however, 
is very laborious, and it will be found that where the 
final error is of the order of that in the above example, 
namely 1 per cent., the results obtained are sufficiently 
accurate for all practical purposes. 

The natural frequency now is 


The values of a, given in Column 3, when plotted 
as shown in Fig. 8, give the “ normal elastic curve ” 
of the system. This curve shows the relative amplitudes 
of the various parts of the system for free vibration, 
and these amplitudes remain in the proportion given, 
whatever the actual value of a may be at mass No. 1. 

The point p (Fig. 8) where the normal elastic curve 
crosses the base line is the node of the system. Thus 
this method gives the positions of the nodes as an 
incidental part of the frequency calculation. 

In Column 8 are given the stresses in each section 
of the shafting due to unit deflection at No. 1 crank, 
the shaft diameter being 8 inches. These values, which 
are filled in after the frequency calculation is completed, 
are obtained from the relationship 

T 
«De 
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where S = torsional stress in lbs. per sq. in. 
T = torque on section in lbs. inches. 
Zr = section modulus. 


3 
= ~ (inch® units), d being the shaft diameter 


in inches. 


The total torque on each consecutive section in ft. lbs. 
for one radian twist at No. 1 crank is given by Column 5, 
so that Column 8 is calculated from 


12 (Value in Column 5) 
57 -3rd/16 


1:06 (Column 5) 
eg 


lbs. per sq. in. per degree. 





Fig. 9.—Variation of Stress in Crankshaft. 


The various stresses are shown diagrammatically in 
Fig. 9. 

If there is a variation in shaft diameter at any section 
of the shaft, say from d to d,, the stress per degree at 


that section will be altered in the ratio (4/ d y’ 
1 


Gorfinkel’s Method.—A modified form of the above 
calculation, in which the work involved is reduced by 
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means of the use of non-dimensional coefficients, has 
been devised by Gorfinkel.* Referring to the last of 
Equations (28) 
@? 
Ant) = a, a R. (1,a, a I.a, — oot it [,a,) 


ie 
— 2 2) 2 
= 4, —w?. R=! ny, 
I l R : 
Lett Se pg, lee as aes 7a here / is th 
etting 0, / L, A, / l, / RY ere lis the 
*““ equivalent length” of shaft used in calculating the 


2 
stiffness R, a, = | and X? = ® ti] R» we have 
1 1 





a 
X*R, iA 2 
1 oy 
h1, a, = “a, — - Nol ,a.a 
at] “1 nt I, R, 5 en 1”), J 
n 
Ayty = %, — Xx? A, % Ont a : : : (30) 


Assuming any value for X%, successive values of « 


are determined, and if X? has been correctly chosen, 
n+1 
200 will equal zero. 


Remembering that X? = w*J,/R, and that the natural 
frequency is given by w/2zx, we have 


R, 


Natural frequency = x 7 VPS - . (81) 
; } 


Again in this case, the calculation is most conveniently 
arranged in tabular form. In what follows, the previous 
example has been worked out by this method. 


“Engineering, 27th December, 1929, p. 827. 
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TABLE II 


A? = 0636 | 











— e a Zea A A*)Zpa 
1. l 1-0000 | 1-0000 0636 
0636 | 9364 | 1 
| l 9364 | 1-0364 1230 
“1230 8134 1 
| I ‘8134 | 22-7498 1748 
: A748] 6386 | 1 
l 6386 | 3-3884 2154 
oust | 4292 | 
1 4232 | 3-8116 -2423 
+2423 -1809 l 
1 1809 | = 3-9925 2511 
‘2511 |—3-9452 | -99 
7. 66-2 | —-0702 0473 
The natural frequency is 
X /R, 1. /0-0636 x 14-1 X 10° 
: SN oan a 56:5 v.p.s. 


It will be noted that the ordinates of the normal 
elastic curve are again given by this method, namely, 
by the third column in the table. 

Approximate Frequency Calculations.—The methods 
of calculation described in the two preceding sections 
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must be used when the natural frequency is required 
to a high degree of accuracy and when, for the 
calculation of stresses and amplitudes, the normal 
elastic curve is required. 

It occasionally happens, however, particularly in the 
preliminary stages of design, that an approximate 
figure for the natural frequency is required. In these 


(a) 


(2) 


Fig. 10.—Six Cylinder Engine and Equivalent Two-Mass System. 


cases it is possible to use abbreviated methods of 
calculation which will give the one node natural 
frequency to within about 5 per cent. 

(a2) The Concentrated Mass Method.—Referring to 
Fig. 10 (a), which is a reproduction of Fig. 8, if the 
whole of masses 1 to 6 is concentrated at the centre 
line of the engine, a simple two mass-one shaft system 
is obtained. This is shown in Fig. 10 (b), The combined 
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moment of inertia now equals 42-66 lbs. ft.* sec.*, the 
equivalent length of shaft is 8-14 ft., and the stiffness 
4-05 x 10° lbs. ft. per radian. 

This gives a natural frequency 


4:05 X 10% x 442-66 


aie BIG ; 
f= 5 400 x 42-66 = wep 


as against 56-5 v.p.s. for the exact calculation, the 
error being 8-7 per cent. 

It will be seen from the above that to group the whole 
of masses 1 to 6 at the engine centre line gives a fre- 
quency which is too low. A much closer approximation 
is obtained by concentrating a certain percentage of 
the masses at the engine centre linc. This percentage, 
which varies with the number of cylinders, is given 
in the following table :— 


Number of cylinders, . =k 2 3 4 6 6 7 8 
Percentage of Total Masses, 100 90 86:5 85 85 85 85 865 

Thus in the example under consideration, the grouped 
mass now becomes 0°85 xX 42-66 lbs. ft.* sec.? and the 
natural frequency 


1 fo x 10° x 486-3 
=a 400 x 36-3 
the error now being 1-77 per cent. 

(b) The Graphical Method.—The degree of accuracy 
of the graphical method depends on three conditions, 
all of which are fortunately fulfilled by the majority of 
present-day diesel engines. These conditions are :— 

(i) The moments of inertia per cylinder are equal. 
4 . moment of inertia of flywheel/ TI, 
(il) “The ratio moment of inertia per cylinder ( - | I, ) 


is not less than 40. 


= 55°5 V.p.s. 
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(iii) The equivalent lengths of shaft between cylinder 
centres are equal. 


Fig. 8 is a typical diesel engine system which fulfils 


KERNS 
ee SE Se 
EEA 


Frequency VPS. 





7 Z 3 4 5 6768691N 


Ratio -—— Tee 


Fig. 11.—Natural Frequency Graph. 


these conditions and Figs. 11 and 12 are the graphs 
from which the natural frequency may be read direct. 

In Fig. 11 the natural frequency f is plotted against 
the ratio R/I x 10° for engines having from 4 to 10 
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cylinders. The ratio l/l, is assumed to be unity and the 
ratio I,/I, is taken as 60, this being a good average value 
for present-day diesel engine practice. (ly 1s the equiv- 
alent length of shaft between the last cylinder and the 
flywheel, J, is the equivalent length between cylinder 
centres.) 

Since the graph is a series of straight lines when 
plotted logarithmically, the natural frequency may be 
expressed as a single equation in terms of the gradient 
and the intercept on the f axis. The gradient is, as 
would be anticipated, 0-5 and the equation therefore 


becomes 
= 103 ae V.p.S 


where U has the following values :— 


10 cylinders. ; . U = 25-45 


9 cylinders ‘ : . 27-9 

8 cylinders. : . 81:05 
7 cylinders . B49 

6 cylinders ; , . 40 

5 cylinders : . 46:9 

4 cylinders : . 56-9 


Fig. 12 is a supplementary graph to take into account 
variations in the ratio lp/l,.. In this graph the multiplying 
factor V has been plotted against the ratio l,p/l, for values 
of the latter from 0-8 to 1-2. Thus the natural frequency 
of any diescl engine complying with conditions (i), 
(ii) and (111) above is given by 


UV |R, 


l= Io® i, V.p.s.  . ; (82) 
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where U has the values already given and V is read 
from Fig. 12 for the particular number of cylinders 
and l,/l, ratio under consideration. 

The Introduction of Gearing.—The introduction of 
gearing into a torsionally elastic system renders the 
problem of torsional vibration more complex than that 
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Fig. 12.—Graph of Correction Factors. 


of the normal straight system. The ratio of the gears 
must be allowed for, whilst if two equal systems 
are mutually geared to a third, two separate and 
distinct series of frequencies are possible. This is 
made clear by examination of the following simple 
systems. 

Fig. 18 (a) shows two masses and shafts mutually 
geared together. The system is assumed to be 
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(a) 








Fig. 13.—Single Geared System. 


constrained except for torsional movement. The 
nomenclature is as follows :— 

I, = moment of inertia of first mass. 

I, = moment of inertia of small gear. 

I, = moment of inertia of large gear. 

I, = moment of inertia of second mass. 

R, = stiffness between first mass and gears. 

R, = stiffness between gears and second mass. 

o = gear ratio, t.e. ratio of smaller to larger gear. 


It can be shown* that the natural frequencies of this 
system are given by the equation 
FF? — F( 9, +j2+0°Js+Ja) +js(0%J3 +e) tJ oje= 9 (38) 
* Bradbury, Proc.J.Mech.H., 1933, Vol. 125, p. 485 et seg. 


NATURAL FREQUENCIES OF MULTI-MASS SYSTEMS 383 


The notation is the same as that for Equation (23) 
and the combined moment of inertia of the gears is 
(I, + o7f,). Thus it is seen that the system may be 
represented by an equation of the form of Equation (23) 
by putting » = 2 and multiplying the moments of 
incrtia and the stiffness beyond the first gear by the 
square of the gear ratio. The equivalent system, in 
which all masses rotate at the same speed, is shown in 
Fig. 13 (0). 

The system shown in Fig. 14 (a) is the same as that 
of Fig. 13 (a), except that the first mass, shaft and gear 
are duplicated. 

It can be shown* that the natural frequencies of this 
system are given by the equation 


(Fo 9,}{F? — F(j, + 292 + 0773 + ja) 
+ 5:(07%3 +34) + 23eJ4} = 90 (38a) 


The notation is again the same as that for Equation (23) 
and the combined moment of inertia of the gears is 
(22, + 0°I;). 

The second bracket of Equation (88a) proves that the 
system of Fig. 14 (a) may be represented by an equation 
of the form of Equation (28) provided that the first 
mass, shaft stiffness and gear mass are doubled, and 
each subsequent mass and shaft stiffness is multiplied 
by the square of the gear ratio. In addition the first 
bracket shows that there is a frequency due to the two 
equal masses swinging one against the other, the points 
of juncture with the gears being the nodes. The 
equivalent systems in which all masses rotate at the 
same speed are shown in Fig. 14 (0) and (ec). 

These two examples serve to illustrate the type of 


* Bradbury, Proc.J.Mech.H., 1983, Vol. 125, p. 485 et seg. 
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vibration to be expected. In the geared diesel engine 
drive, however, the number of masses and shafts is so 
great that it is both quicker and more advantageous 


(a) 





() 





i, 


(c) 


Fig. 14.—Double Geared System. 
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to calculate the frequencies by successive approxima- 
tion. 

This method of calculation is carried out in the 
manner already described in Table I and the calculation 
proceeds in the normal manner until the gears are 
reached. At this point, the angle of twist (Column 8) 
is multiplied by the gear ratio and the total torque 
(Column 5) is divided by the gear ratio. The calculation 
then proceeds as before. The normal elastic curve and 
the total torque on each consecutive section of shaft 
per radian twist at No. 1 mass are again obtained as 
an incidental part of the frequency calculation. 

The following table gives the calculation for the one 
node frequency of the geared marine diese] engine 
system shown in Fig. 15. The specification of the 
system is as follows :— 


Moment of inertia per cylinder 
line. : ; : . = 18-73 lbs. ft.? 


sec.? 
Moment of inertia of flywheel = 388 lbs. ft.? sec.? 
Moment of inertia of first gear. = 2 Ibs. ft.? sec.? 
Moment of inertia of second gear = 42 lbs. ft.? sec.? 


Moment of inertia of propeller 
(including entrained water) . = 174 lbs. ft.? sec.? 
Stiffness between cylinder centres = 18:12 x 108 ft. 
Ibs. per radian 
Stiffness between No. 7 cylinder 
and flywheel : ; 18:12 x 10° ft, 


Ibs. per radian. 


Stiffness between flywheel and 
first gear . . : - = O-815 x 10 ft. 
Ibs. per radian. 


t 
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86 
Stiffness between second gear 
and propeller. ; . = 1:005 x 106 ft. 
Ibs. per radian. 
. = 1/25 


Reduction gear ratio 





Fig. 15.—Seven Cylinder Marine Diesel Engine with Reduction Gearing. 
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The remainder in Column 5 shows an error of 
approximately one per cent., so the value of w = 62:1 
may be taken as correct. 

The natural frequency thus is 


f = wf2x = 9-9 v.p.s. 
= 594 v.p.m. 


In Fig. 15 is plotted the one node normal elastic 
curve as calculated in the above Table, together with 
the two and three node forms which are calculated in 
exactly the same manner. 

Referring back to Fig. 14, it was stated in the case of 
the two node frequency that the two equal masses swing 
against each other, the points of juncture with the 
gears being the nodes. This brings out a very important 
point in connection with twin engines geared to a 
common propeller shaft. 

Provided that the engines are identical in every 
respect and that the cranks are absolutely in phase 
(this is the usual arrangement for synchronised starting 
and reversing) the mode of vibration just described is 
unexcitable since the disturbing forces attempt to 
oscillate the crankshafts in phase and the natural mode 
of vibration is exactly out of phase. 

Thus with this type of installation, torsional vibration 
with an even number of nodes cannot appear in practice, 
and only those with an odd number of nodes need be 
considered. 
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EQUIVALENT DYNAMICAL SYSTEMS 


In all examples discussed so far, what are known as 
‘“‘ ideal’? systems have been considered. That is, the 
masses have been assumed to be uniform circular rings, 
and the shafts weightless and uniform in diameter. In 
actual practice much more complex conditions prevail, 
so that it is necessary, before calculating an actual 
mstallation, to replace it by an equivalent “ideal ” 
system. Fig. 16 shows an actual dynamical system 
comprising a six cylinder diese] engine direct coupled 
to an alternator, together with its equivalent “ ideal ” 
system. 

(a) Equivalent Lengths of Shafts.—By equivalent 
Jength of shaft is meant that length of uniform shaft 
which, when subjected to a given torque, will twist 
through an angle equal to that produced by the same 
torque on the variable section shaft. 

(i) The Crankshaft.—The engine crankshaft is a very 
complicated member, the equivalent length of which 
cannot be accurately calculated from first principles. 
As a result, however, of some very valuable work by 
Major B. C. Carter, an empirical formula* has been 
devised which is found to give very reliable results 
in practice. Referring to Fig. 17, the formula is as 
follows :— 





Equivalent length from A to B in terms of journal 
shaft 


3(D,4 — d,ja ar (D,*—d,*) 
L = (2b + 08h) + 7 pag +e 2 hme (84) 


* Hngineering, 13th July, 1928, p. 36. 
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Corresponding stiffness from A to B 


Gr(D,4 -— d,*) 
82L 





oO 


(35) 


where G = 17 X10® lbs. per sq. ft. for steel. 
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Fig. 16.—Six Cylinder Diesel Engine and Rotor and Equivalent 
Ideal System. 
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This formula was obtained by twisting actual aero, 
automobile and marine engine crankshafts in bearings 
having clearances approximating to those used in 
ordinary working. It will be noted that, in the first 
term, it has been necessary to allow for local yielding 
of the webs at the point of juncture with the journals 
and crankpins by adding the length 0-8h. Also, from 
the second term, it will be scen that the crankpin is 
in pure torsion, 7.é. when the crank is twisted, the 


WT) Ca 





Fig. 17.—Crankshaft (Carter’s Formula). 


journals move out of the axial. Thus when a crankshaft 
is vibrating torsionally, the journals move to and fro 
across the clearance spacc. 

(ii) Variable Diameter Shaft.—Since the angle of 
twist produced in a given uniform shaft by a given 
torque is inversely proportional to the fourth power of 
the diameter, the equivalent length at diameter d, of 
a variable diameter shaft such as that shown in Fig. 18 
is given by 


L=1, +4, (+z) 
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Fig. 18.—Shaft of Varying Section. 


(ill) Zapered Shaft.—For the tapered shaft shown in 
Fig. 19, consider the elementary strip of width 6l,. Let 


the diameter at this point be (d, + fl,) where B is a 
constant. 





Fig. 19.—Tapered Shaft. 


Then the angle of twist in the short length 4, is 


Cél, 
(d, + l,)* 


where C is a constant. 
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Integrating, we have 


Total angle of twist between [ C 
AandB. = 0 (d, + Bl,)* él, 


Ol (att dete + ait) 
~ 3 d,3d 8 





since B = “aa 4, 
Now if £ is the length of an equivalent shaft of 
diameter d,, we have 





CL Cl /d2 + did, + d,? 
d,4~ 3 ( d,3d.3 


whence 
hte 
L=, G + gat gs . (86) 


When making calculations for equivalent length, it 
is essential that the various parts are given their 
correct amount of elasticity. Thus, flanged couplings 
Should be calculated up to the bolt centres only ; 
flanged couplings keyed on to the shaft should be 
calculated on the shaft diameter for half the length 
of the key and on the coupling diameter for the 
other half. Alternator rotors and generator armatures, 
pressed and keyed on to the shaft, should be similarly 
treated. 

(6) Equivalent Masses. 

(i) Moment of Inertia of Flywheel_—Fig. 20 shows a 
section of a typical diesel engine flywheel. It is 
usual, for ease of calculation, to divide the section up 
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eines 





Fig 20 —Flywheel 
into rectangular section rings as shown at 1, 2 and 8. 
Then 
: Te 
Weight of Part I = W, = qia? — d,*)w, X 450 Ibs. 


for cast iron (d,, d, and w, being measured in feet). 


Radius of gyration = r, = fete ft. 


; : W,r/? 
Moment of inertia of Part I = ar lbs. ft.2 sec.? 


Thus the moment of inertia of the whole flywheel is 
zs Wary t+ ae + W,r,? tbs. 


It is usual in these calculations to neglect the flywheel 
effect of the hub. 


{t.? sec.? 


I 
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(ii) Armatures and Rotors.——When calculating the 
moment of inertia of an armature or rotor there is no 
accurate way other than to find the Wr? value of each 
individual part and finally to add them together as in 
the case of the flywheel. Though this method is 
somewhat laborious, extremely accurate results can be 
obtained by careful calculation. 

A method of finding the moment of inertia by 
experiment will be described later. 

(ii) The Engine Parts.—Owing to the fact that the 
reciprocating parts do not have a constant inertia 
value, it is usual to substitute an equivalent revolving 
mass for the combined revolving and reciprocating 
masses of each cylinder. To do this, each part must 
be considered separately and a final summation made 
of the Wr? values. 

The Crank Pin.—Referring to Fig. 17, the weight of 
the crankpin (assumed solid) is 


W, == x Dj x a x 489 Ibs. 


The radius of gyration may, without appreciable error, 
be taken as the crankthrow r. 
For the bore, the weight to be deducted is 


W, =| x d? x (a + 2h) x 489 Ibs. 
The radius of gyration is again r, so that the moment 
of inertia of the crankpin less the bore is 
(W, — W,)r’ 
8g 

The Crankwebs.—The weight of the crankwebs (if 

rectangular) is given by 
W, = 2hws x 489 lbs. 


I= Ibs. ft.? sec.? 
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The radius of gyration is given by 


guile 4 an) (5-») ft. 


(where v is the distance from the centre line of the 
journal to the base of the ae 


and the moment of inertia is ——-— = Abe, {t.? sec.2 


If the crankwebs are not rectangular it will probably 
be more convenient to treat them in the manner given 
below for balance weights. 

The Balance Weights.—These are best calculated by 
a graphical method. Referring to Fig. 21, a typical 
balance weight is shown with its two principal axes 
XX and YY, XX being the crank journal centre. 
The procedure for finding the moment of inertia is :— 

On the right hand half of the balance weight, a 
rectangle OPQR is described. On to the base RQ, 
draw any line OA. Erect a perpendicular to cut the 
edge of the weight at B and C. Draw lines horizontally 
from B and C to cut the line OA at D and E. 

Draw a number of lines similar to OA and follow the 
same procedure. Then by joining up the points D, E, 
etc., what is known as the first derived figure is obtained. 
This is shown in single shading on the right of Fig. 21. 

The next step is to treat this first derived figure in 
exactly the same manner and so obtain the second 
derived figure shown double shaded. 

Then calling the areas of the original, the first derived, 
and the second derived figures’ A,, A, and A, 
respectively, the radius of gyration of the balance 
weight about the axis XX is given by 


A, x (PQ)? 
- REPT oy 
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If the thickness of the weight is h and the material 
cast iron, the weight is given by 


W = 2A,h x 450 Ibs. 


E90 
ZO KS 
EXD KX 


C; 
DORR 


Fig. 21.—Balance Weight. 





and the moment of inertia about the axis XX is 


2 
Wrsx Ibs. ft.2 sec.® 





ly, = 
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It is now necessary to find the moment of inertia 
about the axis YY. To do this, the same procedure is 
followed except that ST (shown on the left of Fig. 21) 
is used as the base. The first and second derived 
figures are then as shown. Using the symbols A,, 
A} and A,! in the same manner as above, the radius 
of gyration about the axis YY is given by 


‘A,1 x (OS)? a 


T = 
YY / A, 
and the moment of inertia by 


Wry? 
——~ bs. ft.2 sec.? 


I yy = 

The rotational moment of inertia about the axis 

perpendicular to the paper through the point O, 1.e. 
the crankshaft centre line, then is 


FT otational = Peg + Ise lbs. {t.? sec.” 


The Connecting Rod is neither wholly revolving nor 
wholly reciprocating, but a very close approximation to 
its true value can be obtained by splitting it up into 
its component revolving and reciprocating parts. This 
is usually done by the two scale method shown in 
Fig. 22. 

The rod is carefully balanced in a horizontal position 
on two scale pans, knife edges being applied to the 
large end and small end centre lines. The weight W,, 
registered by the large end scale is taken as revolving, 
and is added to the crankpin at the latter’s radius of 
gyration. The weight W,, registered by the small end 
scale pan, is taken as reciprocating and is added to the 
piston and gudgeon. 
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The Piston and Gudgeon, which are wholly recipro- 
cating, must, together with the reciprocating component 
of the connecting rod, be replaced by an equivalent 
weight at crankpin radius. Referring to Fig. 28, let 





Fig. 22.—Connecting Rod Weighed by Two Scale Method. 


the total reciprocating weight be W, and the equivalent 
weight W. Then if the crank is revolving with an 
angular velocity w radians per second, the kinetic 
energy of the weight W is 

Wr'w? 


E,, == og ft. lbs. 


w. + 
We oR 






Fig. 23.—Equivalent Revolving Weight. 


Now neglecting connecting rod _ obliquity—which 
entails very little error for connecting rods more than 
four cranks long—the velocity of the weight W, in the 
crank position shown is 


v=or sin wt f.p.s. 


where ¢ is the time in seconds. 
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Thus its instantaneous kinetic energy is 
Wor sin wt)? 
2g 
and the average kinetic energy, taken over one 
revolution, 1s 


ft. lbs. 








l W w?r2 2r 
seas, teas ¢ m2 
i ic a * 2g i » sin wt dt 
— Wor’ 
= ag 
W. 
Thus if E, = E,,, then Toes Therefore one half 


of the reciprocating weight should be added to the 
purely revolving weight at crankpin radius. 
To summarise, the masses concentrated at each 


cylinder centre line are :— 


The crankpin (less bore, if any) at crankpin radius. 

The crankwebs, acting at their radius of gyration. 

The balance weights (if any) acting at their radius of 
gyration. 

The revolving portion of the connecting rod, at 
crankpin radius. 

One half of the reciprocating parts, 7.e. reciprocating 
portion of the connecting rod, piston, rings, gudgeon, 
piston rod and crosshead (if any), at crankpin 
radius. 


(iv) Moment of Inertia by Eaperiment.—The moment 
of inertia of a complicated member such as an armature 
or rotor may be found quite easily by experiment. 
Referring to Fig. 24, the mass is slung by two equal 
length wires from a beam, the axis of the mass being 
vertical, and is given a small initial twist and allowed 
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to oscillate torsionally. The time taken to make one 
complete oscillation is then noted. If this time is 
t seconds 
Ur? 
b= 25 ae: 
gx 
where IJ is the length of each suspension wire in feet. 
z is the distance of each wire from the central axis. 
g is 82-2. 
r is the radius of gyration in feet. 





Fig. 24.—Moment of Inertia by Experiment. 


Thus the radius of gyration is found and, the weight 
having previously been noted, the moment of inertia 


2 
is found from the usual formula J = dan 


8 


52 TORSIONAL VIBRATION IN DIESEL ENGINES 


THE DISTURBING FORCES—CRITICAL 
SPEEDS 


Harmonic Analysis.—On page 10 it was shown that 
resonance conditions can be produced in the torsional 
pendulum (Fig. 1) by oscillating it at the correct 
frequency. In the diesel engine similar conditions 
can prevail, as the torsional elastic system, comprising 
the crankshaft and its associated revolving and 
reciprocating masses, is acted upon by the periodically 
varying cylinder torque. 

This cylinder torque is made up of a large number of 
harmonics of varying amplitude and frequency. Each 
of these harmonics is theoretically capable of producing 
resonance at the appropriate number of revolutions 
per minute, the actual degree of resonance produced 
being dependent on the magnitude of the harmonic 
and the amount of damping present. 

The turning effort curve may be resolved into its 
harmonic components by means of the Fourier analysis. 
Fourier’s theorem states that any periodic function 
such as 

y = 2(8) 
which repeats itself in the intervals 0 to 2x, 2x to 4n, 
etc., may be represented by the series 


y = P(8) = dy + a, sin 6 + a,sin 26 +... a,sin pb 
+ b, cos 6 + b,cos 20 +... 5, cosp0 (87) 
or alternatively by 
y = £(6) = ay + a, sin (6 + 9,) 
+ @, sin (20 + g,) +... (88) 
where a = Va? + b,? 
tan 9; = b,/a, 
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In these expressions, @) represents the mean height 
of the curve, whilst the remaining terms of the series 
are known as the harmonic components. 

To perform the harmonic analysis, the mean height ag 
is first determined by finding the area by planimeter 
and then dividing by the length. To find the harmonics, 
the period is divided into 2h equal parts and ordinates 
erected. 

Then the values of the harmonics are given by* 


1. . SNt 
= pre sin > - 
b, = 524, cos 


where A, is the ordinate at 6 = sx/h. 


It will be seen from the above that the number of 
ordinates must be twice the order of the highest harmonic 
to be extracted. Also, it can be proved that the larger 
the number of ordinates, the more nearly will the 
calculated value of the harmonic approach the true 
value. 

A typical example of harmonic analysis is given in 
Fig. 25 in which the turning effort curve is that of a 
four cycle diesel engine. The fundamental harmonic 
is designated “‘ Order 4” since its frequency is one half 
engine speed, the second harmonic “ Order 1,”’ and 
SO on. 

When carrying out the harmonic analysis of a turning 
effort curve, the inertia forces may be included with 
those of the gas pressure and the resultant diagram 


*For a very complete treatise on Harmonic Analysis see ‘‘Fourter’s 
Theorem and Harmonic Analysis,” by Eagle. . 
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analysed, or the two curves may be analysed separately 
and the appropriate harmonics added vectorially. In 
the example given, inertia has been included. 


; MEAN TURNING EFFORT =43°4 LB. PER SQ. IN. 






a 
© 
@ 
a 
o 


MAX. TURNING EFFORT — LB. PER SQ. IN. 
ORDER OF HARMONIC 


Fig. 25.—Harmonic Analysis of Turning Effort Diagram. 


THE DISTURBING FORCES—-CRITICAL SPEEDS 55 


To evaluate the inertia harmonics separately, it is 
not necessary to analyse the inertia curve graphically 
as in the case of the gas pressure. All that is necessary 
is to set down the expression for the turning effort due 
to inertia, and expand it in a series. This series will 
be found to take the form 


J 
Ca sin 0 + ¢, sin 20 


& 
+c,sin 86 + ...c,sinp0) . (89) 





y= 


i.e. there are no cosine components and no half order 
components. 
In the above expression, 


W, = weight of reciprocating parts in lbs. 
g = 82-2 
r = crank radius in feet. 
w = angular velocity of crankshaft in radians 
per second. 


Tables IV and V give the harmonic components for 
gas pressure and inertia respectively when the curves 
are analysed separately. The gas pressure harmonics 
are for a mean indicated pressure of 90 lbs. per sq. in. 
Both sets of harmonics are for a crank/connecting rod 


ratio of 1/4°5. 
The harmonic torques due to gas pressure then 


are :— 
Sine Component = H, =a, X nd?/4 xX r ft. Ibs. 
Cosine Component = H, = b, X nd*/4 X r ft. lbs. 
Resultant Component = H, = 2, X nd?/4 x r ft. Ibs. 


where d is the cylinder diameter in inches and r is the 
crank radius in feet. 
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56 
TABLE IV 
Fiaeawnier | ~ Sine ft Gone Revultane: 
Order n Term a, Term 6, Term z, 
4 222 —23 32 
35 


32-5 —13 


1 
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The harmonic torque due to inertia is :— 
Sine Component = H, = c, X W.w?/g x r? ft. Ibs. 


The resultant harmonic torque due to the combination 
of the gas pressure and inertia harmonics is given by 


H,, = V(H, + H.)? + H,? ft. Ibs... (40) 





TABLE V 

Harmonic Order n Sine Term c, 

1 + -056 

2 — +500 

3 —-170 

4 —-013 

5 + -002 

6 +-0001 

7 — -00002 

8 0 


Critical Speeds.—When the crankshaft is revolving at 
such a speed that one of these harmonics synchronises 
with the natural frequency, resonance occurs, and it 
is said to be at a critical speed of that particular order. 
Thus in a four cycle engine it is possible to have critical 
speeds of every integral order and half order and in 
two cycle engines, by similar reasoning, critical speeds 
of every integral order, the order number being deter- 
mined by the number of vibrations occurring per 
revolution of the crankshaft. 
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Now it will at once be apparent that the magnitude 
of the vibration in a multi-cylinder engine differs widely 
fur the various critical speeds, even when allowance is 
made for the difference in value of the harmonic forcing 
torques. This is due to the fact that, for some orders, 
the harmonics assist one another in producing vibration, 
whilst for others, the reverse process occurs. It therefore 
becomes necessary, not only to classify the critical 
speeds, but to obtain some idea of their relative capacity 
for producing vibration. 

Let an nth order harmonic torque H, be applied at 
the rth cylinder where the ordinate on the normal 
elastic curve has the value a,. 

The work done, 2.e. the kinetic energy, at this point 
is proportional to 

Ha, 


Then for the entire group of cylinders, assuming the 
torque H, to be applied at each cylinder centre line, 
the total kinetic energy is given by 


E =a constant X H,2a, . ~ (41) 


where 2a, represents the vector sum of the normal 
elastic curve ordinates at which the torque is applied. 
We thus reach the important conclusion that the 
magnitude of the vibration in a critical speed is directly 
dependent on the value of H,,2a,, that is, on the harmonic 
forcing torque, on the shape of the normal elastic 
curve, and on the crank sequence and firing order. 
The value of 2a is most easily evaluated by means 
of phase and vector diagrams, a typical example of 
which is shown in Fig. 26. The system is that of a 
six cylinder four stroke engine with a firing order of 
1, 5, 8, 6, 2, 4. The normal elastic curve is as shown. 


and 2% 


WIN 





Zand 2 


Fig. 26.—Phase and Vector Diagrams, Normal Six Cylinder Engine. 
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The explanation of the diagrams is as follows. For 
each order of critical speed, a phase diagram is con- 
structed, 360° representing one complete vibration. 
For the half order critical speed, the cylinders come 
into action in the phase diagram at equal intervals and 
in their correct firing order as shown at (a). For the 
first order critical speed, since one vibration is completed 
in half the time, the angles will be doubled as shown 
at (b). For orders 14, 2, etc., the angles for order 4 
are multiplied by 3, 4, etc., up to and including order 8. 
After order 3, the diagrams repeat, i.e. orders 4, 34, 
64, etc., orders 1, 4, 7, etc., have the same phase 
diagrams respectively. 

The vector diagrams are shown beneath the phase 
diagrams in Fig. 26. The vectors are drawn parallel 
to their respective cylinders in the phase diagrams and 
of lengths equal to their respective ordinates in the 
normal elastic curve. The closure of each vector 
diagram gives the value of 2a, for that particular 
order or orders. 

It will be noticed that the phase diagrams for orders 
4 and 2} and for orders 1 and 2, respectively, are 
mirror images. Only one vector diagram is therefore 
needed for each pair. 

It will now be seen that for critical speeds of orders 3, 
6, 9, 12, etc., the harmonic torques are all in phase 
and tend to oscillate the shaft in its normal mode of 
vibration. These critical speeds are the most powerful 
and are known as “major” critical speeds, all others 
being termed “‘ minors.”” For the type of system under 
consideration, which incidentally is the type most 
frequently encountered in diesel engineering practice, 
major critical speeds may be defined as those whose 
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order numbers are equal to or multiples of the number 
of power strokes per revolution. 

It should be clearly understood that the above 
definition of major critical speeds does not apply to 
all systems. Whether a critical speed is major or 
minor depends entirely on the shape of the normal 
elastic curve. This will be made clear by reference 
to Fig. 27, which shows the normal elastic curve, phase 
and vector diagrams for a system which differs from 
that of Fig. 26 only in that a flywheel has been placed 
either end. The node, in this case, lies in the centre 
of the engine crankshaft and the vector diagrams show 
that for such conditions, orders 1}, 44, 74, ete., are the 
major critical speeds, all others being minors. 

The Effect of Firing Order.—In the previous section 
it was shown that the firing order determines the 
sequence of the vectors in the phase and vector diagrams. 
Hence the effect of the harmonic forcing torques, as 
represented by the closure of the vector diagram, can 
be reduced or increased for certain critical speeds by 
a suitable modification to the crank sequence and 
firing order. 

This method of modifying the magnitude of critical 
speeds is, however, very limited in its application. 
Crank sequence can very rarely be altered without 
seriously upsetting engine balance and only in certain 
cases can the firing order be altered without altering 
the crank sequence. The only cases where the method 
is extensively used are those of six and eight cylinder 
four stroke engines, in which several firing orders are 
available for the same crank sequence. 

The effect of altering the firing order is_ best 
demonstrated by an actual example. The equivalent 
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Fig. 27._Phase and Vector Diagrams. 


THE DISTURBING FORCES-— CRITICAL SPEEDS 638 


dynamical system is shown in Fig. 28, together with the 
normal elastic curve and phase and vector diagrams 
for firing orders 1, 5, 3, 6, 2, 4, and 1, 2, 4, 6, 5, 8, 
respectively. The engine, which has a one node natural 
frequency of 9,050 v.p.m., is normally expected to run 
at 1,200 r.p.m. 

For the above natural frequency, the critical speeds 
in the neighbourhood of the running speed (obtained 
by dividing the natural frequency by the order number) 
are :— 


Critical Speed Order 
1510 r.p.m._. : ‘ : . 6 
92, lw wt ti ttttCS 
1290 ___—sié—“»; ; : : : ae 2 
1205 , ©. . . «.  . % 
1130 so, ; : : é . 8 
1065 4 ©. . «.  «. . 8h 
1005 _ sé, : ‘ é ; ~ 2 


We thus have 6th and 9th order major criticals at 
1510 and 1005 r.p.m. respectively, whilst there is a 
74th order minor critical practically on the running 
speed. 

It has previously been shown that the inagnitude of 
the vibration in a critical speed is proportional to 


H,, a, 


where H,, is the harmonic forcing torque and 2a, is 
the closure of the vector diagram. The system of 
Fig. 28 has therefore been analysed on this basis, 
values of H, 2a, having been found for each particular 
order. 
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Fig. 28.—Phase and Veotor Diagrams for Different Firing Orders. 
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The work is best arranged in tabular form as shown 
in Table VI. This table gives the values of H, 2a, 
for the two firing orders in question. 

Thus, from the fifth column, for a firing order 1, 
5, 3, 6, 2, 4, the 74th order minor critical, which lies 
on the running speed, is almost as important as the 9th 
order major. Turning now to the seventh column, by 


TABLE VI 


Firing Order Firing Order 
1, 5, 3, 6, 2,4 | J, 2. 4, 6, 5, 3 

















| 

| Critica] Npeed Order} H, 

| a05: | =a, |. Say are 
1510 r.p.m. 4 6 | 4:5 | 3 9065 | 17-5 3°9065 | 17-5 pe 
1392_—Ci,, ; 63 | 3-49} -5630] 1-96 | 1-2450)} 4-34 | 
1290 _—i,,, : : 7 | 2-80} -2370 -663 | -2370 -663 

| 1205 ss, d 74 | 2-23 | 11-5845 | 3-54 1095 245 
1130 =z, , ‘ 8 | 1:78} +2370 ‘422 | +2370 422 
1065 _—sé=,, : : 8% | 1-45] +5630 816 | 11-2450 1-80 | 
1005 sé, : ; 9 | 1-20) 3-90865| 4-7 3:9065 } 4-7 | 


A Gr CERNE RGR AE AF AOS — ee ~ err me = 


altering the firing order to 1, 2, 4, 6, 5, 8, the 7$th 
order critical almost disappears and the 64th and 8$th 
become of importance. For a running speed of 
1200 r.p.m., therefore, the firing order 1, 2, 4, 6, 5, 8 
is preferable. 

As would be expected from the phase and vector 
diagrams, the major critical speeds are unaffected by 
firing order. 
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VIBRATION AMPLITUDES IN CRITICAL 
SPEEDS 


Energy Balance.—-On page 5 it was assumed, for mathe- 
matical convenience, that damping was proportional 
to the velocity of the vibrations. By this means the 
equation of motion was made linear and expressions 
were obtained for the equilibrium amplitude and the 
dynamic magnifier and the resonance curves of Fig. 3 
were constructed. 

In the case of a multi-mass system, of which the 
diesel engine is a particular example, the inclusion in 
the equations of motion of a damping factor, even in 
its simplest linear form, leads to appalling mathematical 
complications, and the labour involved in obtaining a 
solution is not justified by the accuracy of the results 
obtained. Fortunately, however, the condition of 
resonance permits a different treatment of the problem 
and allows non-linear damping conditions to be taken 
into account without much difficulty. 

The principle of the method used is that during 
steady conditions of resonance, the energy input per 
cycle due to the forcing torques must be equal to that 
absorbed by damping. 

Damping may arise from a number of sources and 
may be either internal or external, that is, in the engine 
itself or in the driven machinery. In the case of external 
damping, the most usual sources are hysteresis losses in 
the driven shafting, and losses due to torque variation 
in a driven member such as a marine propeller, electric 
generator, etc. In the case of internal damping, the 
losses are divided between hysteresis and oil film 





VIBRATION AMPLITUDES IN CRITICAL SPEEDS 67 


dissipation at the main bearings, the latter probably 
being the more important. 

Work Done per Vibration at the Cylinder Centre Line.— 
Let a harmonic torque of the form 


T =H sin ot 


be applied at the cylinder centre line. Then the motion 
produced will be of the form 


a =a sin (wt — e) 


Fig. 29 shows the torque and the motion plotted 
vectorially, a lagging behind H by the angle «. 


Fig. 29.—Vector Diagram of Torque and Motion. 


Let the work done per vibration be 6w. Then 
bw = [Td 
| 


a [ow e 
=[" Haw sin wt cos (wt — e)dt 
0 
= nwtHa sine. ; : ‘ . (42) 


During conditions of resonance, ¢ = 7/2 and therefore 


sin ¢ = 1 so that 
dw = nila é , . (48) 
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Work Absorbed by Damping at the Cylinder Centre 
Line (Damping Pioportional to Velocity).—_ Again let 
a be the maximum amplitude at the cylinder centre 
line. dx 
Let the damping torque be Ma where yw 1s the 
coefficient of friction, and « a sin oat. 

Then for the work absorbed per vibration 


2x da 
Ow, =f. oy da 
2 1a 
= / Mara? cos? at dt 
0 


= nao. ; : : . (44) 


Extension to Multi-Crank Engine.—Coming now to 
a multi-crank engine, the total work done at resonance 
by the nth order harmonic torque will be, as shown 
on page 58, 

6w=tH,2a . 3 - (45) 


or, if 6 be the measured amplitude at No. 1 cylinder 


centre line, where the ordinate of the normal elastic 


curve is unity 
6w=7rH,62a_ . : . (46) 


Similarly, the work absorbed by damping, assumed 
to be equally distributed between the cylinders, is 
given by 

6w, = Tywh*Za* . . (47) 


Now at resonance dw equals éw,, therefore 


tH, 62a = ruwb?Za* 
whence 
_ H, =a 
be Sa? 


. . (48) 
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This expression enables the damping factor for a 
given engine to be calculated in terms of the measured 
ainplitude at No. 1 cylinder and vice versa. 

Referring back to page 9, it was shown that for 
resonance conditions in the torsional pendulum, the 
dynamic magnifier is 1/K, which equals I@/u. This 
same conception can be uscfully applied to the multi- 
crank engine having equal moments of inertia per line, 
as follows. 

Dividing both sides of Equation (48) by J@ and 
inverting, we have 

i lw = Iw?6Xa? 

Dynamic Magnifier - ic Se : . (49) 

Damping Within the Engine---Shannon’s Method.—In 
a paper published in the Journal of the Royal Technical 
College, Glasgow, 1933, and later read in an amplified 
form before the Institution of Mechanical Engincers,* 
Mr. J. F. Shannon published the results of his researches 
on torsional vibration, and he came to the conclusion 
that the vibrational energy at resonance was dissipated 
mainly by the induced actions on the oil films at the 
big end and main bearings. Losses due to hysteresis 
were found to be such a small percentage of the total 
losses that they could be neglected. Therefore by 
assuming that the overall damping was viscous in 
character, the dynamic magnifierf was correlated to 
the vibration form of the engine and a relationship 
was obtained whereby the amplitude in a critical speed 
could be reasonably determined. 


*** Engine Damping Factors in Torsional Oscillation and the Effect of 
Vibration Form,” by J. F. Shannon, J/.Roy.Tech.Coll.,Glasgow, January, 


1933. 
¢ Called by Shannon “ the Engine Damping Coefficient.” 
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It was found that for different engine systems, there 
was a wide variation in the value of the dynamic 
magnifier. The latter was, however, found to be 
dependent on the shape of the normal clastic curve, 
so that by plotting Iw/y against Ya/n, where n is the 
number of cylinders, a reasonable classification of a 
number of test results was obtained. From these 


0 01 02 63 04 05 06 OF O08 09 


Average Relative Amplitude Za/n 
Fig. 30.—Shannon’s Line. 


results, the curve shown in Fig. 80 was put forward 
as a guide in design. 


The method of using the curve is best shown by the 


following worked example of a six cylinder diesel 
engine. The specification is as follows :— 


Moment of inertia per cylinder 
line 


; : : 1-2196 lbs. ft.? sec.® 
Moment of inertia of flywheel 


22°6 lbs. ft.? sec.? 





1 
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Equivalent length between 
cylinder centres g . = 1°58 ft. of 43” dia. 
shaft. 
Equivalent length between 
No. 6 cylinder and flywheel = 1-657 ft. of 43” dia. 
shaft. 


The equivalent system as set out for checking is as 


shown in Fig. 31. 


= 
oe 


2:592x10|2-592 x10°|2-592x10 | 2-592%105| 2-392x105| 2-470x10° 


= 


1-2196 12196 V2196 — 2196 12796 =7-2796 


Leb 
Fig. 31.—Six Cylinder Diesel Engine. 


Using the Successive Approximation method of 
calculation, the table is as follows (see page 72) :— 
The natural frequency 


_ 893 62-6 
gh V.p.s. 


and the normal elastic curve is as shown in Fig. 81. 
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It is now required to find the amplitude in torsional 
vibration due to the sixth order critical speed, which 
occurs at 


62°6 x 60 
——__——-— == 626 I.p.m. 
6 
The ‘‘ form factor ”’ 
3°7364 
Ya/n =" "= = 0-6227 


From Fig. 30, the dynamic magnifier Iw/u corresponding 
to this form factor is 28. 





Hence 
Iw Iw?6Sa? 
bo Oe H,da 
so that 
28H, a 
Iw? za? 


The engine in question has a bore of 63 ins. and a 
stroke of 9 ins., so that 
4°50 x 31°92 X 4°5 
ae 
28 xX 54 X 3°73864 
~ 1-2196 x 154500 x 2-9555 
== 0:0102 radian = 0-584 degree. 


= 54 ft. Ibs. 


el 


From the last column of Table VII the corresponding 
stress at the node will be 0-584 x 3-15 = 1-84 tons per 
Sq. in. 

The actual measured amplitude at No. 1 cylinder 
was 0:-437°, the corresponding stress at the node being 
1:87 tons per sq. in. 

Hysteresis Damping.—The shafting material of a 
diesel engine and its associated drives, when subjected 
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to repeated cycles of stress as in torsional vibration, 
does not obey Hooke’s law exactly. The stress-strain 
relationship of steel under these conditions is generally 


Stress 


otrain 


Fig. 32.—Hysteresis Loup. 


as shown in Fig. 32, the area of which gives the amount 
of energy dissipated per cycle per cubic inch of material. 

When the stress is low, the area of this “ hysteresis 
loop ” is so small as to be negligible. As the stress 


fnergy Dissipated 


dtress Kange 
Fig. 33.—Typical “‘ Energy-Stress ’’ Curve. 


range is increased, however, the area _ increases 
enormously, so that the curve of energy dissipated, 
plotted against stress range, is as shown in Fig. 388. 
A number of authorities have investigated the question 
of hysteresis loss in crankshaft steels. The most recent 
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and probably the most comprehensive series of 
experiments were conducted by Dr. Dorey, who 
presented his results in a paper to the Institution of 
Mechanical Engineers in November, 1932.* Dr. Dorey 
found that all the stecls tested showed a quasi-critical 
point below which the rate of dissipation of energy 
varied approximately as the third power of the stress 
range and above which the rate varied approximately 
as the tenth power of the stress range for primary 
loading and as the sixth power of the stress range for 
secondary and subsequent loadings. 

The following expression for hysteresis energy per 
unit volume per cycle was evolved :— 


S\yr 
ow — 2(5 ) . oy (50) 


where dw = hystcresis energy per unit volume per cycle 
(inch Ibs. per cubic inch). 
Z == hysteresis constant. 
S' = half the stress range for alternating stress. 
S, = half the critical stress range. 
y = hysteresis cxponent. 


For a uniform hollow circular shaft, of external and 
internal radii 7, in. and r, in., respectively, and length 
1 in., the hysteresis loss per cycle is 

tom gg tras 
(2 +) ry” S, 
which becomes 


tn. (8)) 
ho = ep 2-5 Lo. (81) 


for a solid shaft of radius 1. 


* “Elastic Hysteresis in Crankshaft Steels,’’ Dorey, Proc.J.Mech.E., 
1932, Vol. 123, p. 479. 
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For a normal crankshaft material of 0-21 per cent. 
carbon steel having a maximuin tensile stress of 30 tons 
per sq. inch, 


(a) For primary loading 


ZL 0°18, 
S 6 tons per sq. inch. 


y 2-8 for values of S below S.. 
10-0 for values of S above S. 


(b) For secondary and subsequent loadings above the 
critical stress range 


Z — 0°71 

S 6 tons per sq. inch. 

y — 3-2 for values of S below S.. 
6-6 for values of S above S'. 


As the endurance linut for crankshaft steel subjected 
to repeated alternating stress is in the neighbourhood 
of the value of S,, category (a) above may be said 
to represent hysteresis conditions for torsional vibration 
problems in diesel engines. 

Thus for a hollow circular shaft 


O “18% re rae ro 8 5) oe 8 
= —,- 2 .(-) .t. (52 
Ow oa? 28 ( l (52) 


and for a solid circular shaft 


. S 2.8 
bij oe ge, (3) -Linch Ibs. . (88) 


2°4 





Ker Wilson’s Method.—In an article in Engineering 
of February 16, 1984, W. Ker Wilson obtained a 
simplified expression, based on Dorey’s researches, 
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which gives the stress in a critical speed when hysteresis 
is assumed to be the principal damping influence. 
The expression is as follows :— 


S — 602 ‘A/S, Ibs. per sq. in. - (54) 


for 30 ton plain carbon stecl not previously subjected 
to cycles of stress above the critical range and for 
values of S below S,, where 


S is the maximum vibration stress. 


6 


So is the “ equilibrium ”’ stress. 


The value of S, is obtained from a formula for the 
“equilibrium amplitude ” a9, which, converted to the 
symbols used in this book and applied to the type of 
engine under consideration, is 


57°3H,, a 


ay = w?S(Ia?) degrees ° ° . (55) 
where Oey Ss a ee es Ag 
la’: — f,a,? + I,a,4 + oo ec we I,a,?. 


Again using the example of Fig. 31, 


a7°3 X SL X 3°T364 


ba fae eee Bases se. ti ro 
154500 » 4-52 te 


ao = 


Whence the equilibrium stress 
S, = 0:0165 < 3-15 = 0-052 tons per sq. in. 
= 116-5 lbs. per sq. in. 
The maximum vibration stress 


S = 602 '\/116-5 = 8464 lbs. per sq. in. 
= 8-78 tons per sq. in. 
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Summary 
Actual Shannon’s Ker Wilson 
Measurement. Method. (Hysteresis). 
Stress (tons/sq. in.) . 1-37 1°84 3°78 
Correction Factor er 0-75 0°363 





General Comments.—From the results summarised 
above, the following conclusions may be drawn :— 


(a) The “Shannon” method of calculation gives a 
reasonably close approximation to the actual 
measured stress. 


(b) Hysteresis accounts for about one-third of the 
total damping present. 


It must not be assumed from the conclusions reached 
above that hysteresis plays no part in engine damping, 
but numerous investigators* have found that it is 
necessary to multiply hysteresis losses by a factor 
ranging between 2 and 3 to make their calculations 
agree with the results obtained from torsiograph records. 
The Shannon method, whilst apparently ignoring 
hysteresis, actually takes it into account in the general 
over-all damping factor, the greater portion of which 
is found to be due to the viscous losses already 
mentioned. The author’s experience is that the above 
conclusions are applicable to the majority of diesel 
engines. 

External Damping—Hysteresis.—The energy absorbed 
by hysteresis in the driven shafting is calculated by 
means of Equation (52) or (58) by inserting the 
appropriate value of S from Column 9, Table VII. 
This gives the energy loss per degree twist at No. 1 


* See discussion on Dorey’s Paper, loc. cit. 
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crank. If the shaft is of variable section, care must be 
taken to assess the stress and consequent loss in each 
individual section. 

Damping Effect of a Marine Propeller.—J. ¥. Shannon* 
quotes the following formula for calculating ie damping 
factor of a marine propeller 

34-47 


a . . . (56) 


for normal propellers and slips where 


T = propeller torque in ft. Ibs. 
N = r.p.m. 
== damping factor in ft. lb. sec. units. 


Damping Effect of an Electric Generator.—Lewist 
gives the following approximate formula for calculating 
the damping factor of an electric generator. 


422-7E?1, 


MoE B)Ne °° + (87) 


where E = total E.M.F. of generator (terminal voltage 
plus internal drop). 
KE, = back E.M.F. of driven machinery (not an 
external resistance). 
I, = current. 
N ==Y.p.m. 
f = damping factor in ft. Ib. sec. units. 


The method of incorporating these various damping 
factors is again best shown by means of a worked 
example. Stated briefly the method consists of finding 
the energy dissipated in the engine and in the driven 


* Toc. cit. 
ft Lewis, Trans. Amer, Soc. Naval Arch.and Mar. Eng., New York, 1925. 
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unit each in terms of the amplitude 6 at the forward 
end of the engine. The sum of these is equated to the 
input energy from the harmonic forces—also in terms 
of 6—and the resulting equation is solved for 0. A 
complete example of a diesel electric sect is worked out 


in the foHowing section. 


EXAMPLE—DIESEL-ELECT RIC SET 


Fig. 34.—Diesel-Electric Set with Normal Elastic Curve. 


The equivalent system as sct out for checking is as 
shown in Fig. 34. The specification is as follows :— 


Bore and Stroke 

Moment of inertia per cylinder 
line ; 

Moment of inertia of flywhcel 

Moment of inertia of gencrator 

Stiffness between cylinder 
centres 


come 


I 


1l in. X 144 1n. 


6:63 lbs. ft.? sec.? 
376 lbs. ft.? sec.? 
298 lbs. ft.? sec.? 


14-1 xX 108 ft. lbs. 
per radian. 


VIBRATION AMPLITUDES IN CRITICAL SPEEDS 81 


Stiffness between No. 6 
cylinder and flywheel . == 14:25 x 108 ft. Ibs. 
per radian. 
Stiffness of gencrator shaft . — 18-0 x 10® ft. Ibs. 
per radian. 


Using the Successive Approximation method of 
_calculation, the table is as follows (page 82) :— 


269 
Natural frequency = on 7 483 v.p.s. 


There is thus a Gth order critical at 


43 X 60 
pare pe oe 430 r.p.m. 


The stresses in the final column are calculated for 
an 8 in. journal shaft. 

The stress per degree in the 7% in. dia. generator 
shaft will be 


8 \3 
11:5 X (5) == 14-6 tons per sq. in. 


It is required to find the amplitude due to this 6th 
order critical speed. 

Let 6 radians be the amplitude at No. 1 cpanic 
Then, from Equation (46) :— 


The input work per vibration 
= nH,0=a 
= X 4:5 X 95:0388 < 7°25 X 4°87350 
== 475000 inch lbs. 


Damping is divided between three sources, namely, 
(a) Overall damping in the engine. 
(b) Hysteresis damping in the generator shaft. 


(c) Damping in the generator. 
6 
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These dampings will be dealt with in turn. 


(a) Engine Damping (Shannon Method). 


Za = 4°8735. 2a? S24 1225: 
=a 
— _= 0-8122. 
6 0-812 
Whence from Fig. 30, Iw/u = 19-2. 
6°63 x 269 : ; 
Therefore 4 = oo = 92-8 ft. lbs. per radian per 


sec. 
= 1113 inch Ibs. per radian per sec. 
From Equation (47), the work absorbed by damping 
in the engine 
— trum? Xa? 
=m X 1113 X 269 K 4:°12256 
== 3,880,0006? inch Ibs. 


(b) Hysteresis Damping in the Generator Shaft.—This 
Shaft is 7% ins. uniform diameter and is 17 ins, long. 
The stress per degree twist at No. 1 cylinder is 14-6 tons 
per sq. in., whence the stress for @ radians at No. 1 
cylinder 1s 

14-6 X 57°83 X 0 = 8376 tons per sq. in. 


Then from Equation (53) 


: 0-18 S 2-8 
Hysteresis loss = pee ee =) ol 
2-4 ; 


_O18m oe YT ‘3 
oa * ( 6 . 
= 54,500,00002'8 inch Ibs. 





(c) Generator Damping. 
a Fee, 
e  On(E — E,)N? 
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The set was loaded by means of resistance plates, 
hence £, = 0. 
The load conditions were :— 
E = 440 volts, I, = 500 amps., N = 430 
422°7 X 4407 x 500 y 
2m xX 440 < 43802 
= 961 inch lbs. per radian per sec. 
Work absorbed (from Equation (44)) 
== 1 X 961 X 269 x (0°5780)? 
== 272,0006? inch Ibs. 
Equating input to work absorbed, we now have the 
following :— 
475000 = 4,152,00002 + 54,500,00002-5 
whence 54500018 +. 41520 — 47-5 





whence yp = 12 


The solution of this equation, obtained graphically, is 
6 = 0:0088 radian 
= 0-504° 
The actual measured amplitude on test was 0°45°. 


It is of interest to construct a balance sheet of the 
damping. This will read as follows :— 


erp iec — a ee a ee _—— —_—— — -_—~ a enero eee 


Energy Absorbed per 


Location of Damping. Vibration (inch Ibs.). Percentage. 
Engine . : : 300-7 71-95 
Generator Shaft ‘ 96°2 23-0 
Generator ; ‘ 21:1 5:05 


Pence aaTS 


418-0 100-00 








GENERAL NOTES ON CRITICAL SPEEDS 


(i) A cRITICAL speed of appreciable magnitude usually, 
but not invariably, advertises its presence by a rumbling 
noise. This noise is caused by the various moving 
parts—crank journals, pins, gudgeons and_ pistons— 
hammering from side to side in their respective 
clearances, and by backlash in gear teeth. Whether 
an unusual rumbling noise is, or is not, due to a critical 
speed, can be determined by running the engine over- 
speed, cutting out the fuel pumps and allowing the 
engine to “run down.” If the engine runs quietly 
over the greater part of the range, and the noise occurs 
at clearly defined points, then the trouble is almost 
certainly due to resonant torsional vibration. 

(ii) The noise produced in a critical speed depends 
on the temperature of the lubricating oil. An engine 
may run quietly over its whole range for a short time 
after starting up, but when the lubricating oil has 
warmed up and become much Icss viscous, the noise 
can appcar. 

(ili) Noise is no criterion of the stress set up in the 
shafting. An engine with the flywheel located at the 
camshaft end is invariably less noisy in a critical speed 
than the same engine with the flywheel at the opposite 
end to the camshaft drive, yet the stress may be cqual 
in each case. While engines of the former class may run 
quite safely with torsional stresses of +2-5 to +3-0 tons 
per sq. in., engines of the latter class must have these 
stresscs reduced to +1:0 ton per sq. in. on account 
of noise. 

(iv) The amount of noise produced in a critical speed 
is also dependent on the shape of the normal elastic 
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curve. Fig. 35 shows two typical normal clastic curves 
fur diesel generator sets. In Example (a) the motion 
of the crankshaft is great in comparison with the 
generator end of the system, whilst the reverse is the 
ease in Example (b); (a) will therefore be noisy on 
account of the movement of the engine parts in their 
respective clearances, while (b) will be comparatively 





Fig. 35.—Types of One-noded Normal Elastic Curves. + 


quiet as the greater part of the movement takes place 
in the generator where purely torsional movement has 
no effect on bearing clearances. 

In two particular cases which were recently 
investigated, a dicsel-electric sect of type (a) was ex- 
cessively noisy in a sixth order critical speed when 
the measured stress was -+1°5 tons per sq. in. The 
second, of type (b), had a stress of +6-°3 tons per sq. in. 
in the generator shaft, and the critical speed was almost 


inaudible. 
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Fig. 36.—Effect of Critical Speed on Speed Variation. 
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(v) The effect of a critical speed on governing is of 
great importance. An engine, on full load, running at 
a critical speed, will run fairly steadily, yet, immediately 
the load is dropped off, the governor will surge violently. 
This fact is explained as follows. A critical speed 
requires energy to maintain it, and this encrgy is taken 
from the normal developed output of the engine. In 
addition, the energy reaches a peak at the critical 
speed, falling away on either side. The governor is 
therefore faced with the difficulty of controlling the 
engine on a load which falls away whether the speed 
rises or falls. It is incapable of performing this duty, 
and a surge is the result. The trouble is most noticeable 
at light load, as the “critical speed energy ” is an 
appreciable proportion of the normal developed energy. 
At full load it is negligible. 

Fig. 36 shows a series of tachograph records which 
were taken at the same time as the torsiograph test on 
the engine quoted in the example on page 70. The 
governor surge at light load, in the sixth order critical 
at 626 r.p.m., is clearly marked. 


VIBRATION MEASURING INSTRUMENTS 


THE measurement of torsional vibration is an essential 
part of the complete investigation of vibration problems 
in diesel engines. Lack of information concerning the 
various damping factors renders impossible the pre- 
diction of vibration amplitudes from first principles, 
and it is only by “‘ knowing the answer,” in the form 
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of direct measurement, that to-day’s knowledge of 
damping is greater than it was some years ago. 

The instrument used is known as a torsiograph or 
torsional vibrograph and is actually a type of seismograph 
arranged to operate torsionally. The general principle 
is shown in Fig. 37. A light pulley a, is driven from 
the engine crankshaft in such a way that it faithfully 
follows the periodic angular displacements of the 
latter. The flywheel 6b, mounted on ball bearings ‘and 





Fig. 37.—Essential Parts of the Torsiograph. 


connected to a by a light spring or springs, has a sensibly 
constant angular velocity. By a system of links the 
relative motion between a and b is recorded, the record 
so obtained being a representation of the torsional 
vibration together with any other uneven motion that 
may be present. 

There are a number of types of torsiograph in 
existence, but it is proposed to describe here only those 
two types which are commercially in production and 
which are in general use in this country, viz. the Geiger 
and the Cambridge. 
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Fig 38.—Geiger Torsiograph (Operating Side). 





Fig. 39.—Geiger Torsiograph (Pulley Side). 
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The Geiger Torsiograph.—This instrument is illus- 
trated in Figs. 38 and 89, while Fig. 40 is a sectional 
view of a part of the Torsiograph showmg the principle 
of operation. It consists of a very light Aluminium 
Belt Pulley, inside which, and concentric thereto, runs 
a heavy Flywheel. Both Pulley and Flywheel run 
independently on Precision Ball Bearings but are 
connected to each other by a carefully calibrated Volute 
Spring. 


Alurnnum Polley 






B all o 
Beanng if 


- 
eh Hollow 
Spindle 


Fly Wheel 


Fig. 40.—Geige: Torsiograph (Scc tion) 


When the Torsiograph is in use, the pulley is belt 
driven from the shaft, the characteristics of which have 
to be revealed. This pulley thus takes up all the 
torsional vibrations and irregularities of the shaft, while 
the flywheel tends to run at a uniform speed. It will 
thus be seen that there are angular movements relative 
to pulley and flywheel and these are transmitted by 
levers to the tracing pen which marks the diagram on 
a running strip of paper. 
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The tracing lever is provided with different alternative 
fulerums so that the diagram may be obtained in 
different aniplifications according to speed and frequency. 


se vp.m | Zrevs | 


Fig. 41.—Geigcr Torsiogram. 


The paper band is moved either by gearing from the 
pulley or by clockwork. 





Fig. 42.—Cambridge Torsiograph. 


For correctly timing the torsiograms, there is a 
spring, calibrated to make a certain number of 
vibrations (usually 1500 per minute) which are marked 
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on the paper through a second pen. A third pen, 
operated by a contact breaker, cither from the injector 
or the exhaust valve, indicates the revolutions per 
minute and the phase of the vibrations. 

A typical Geiger torsiogram, on which the tracings 
from the three pens are recorded, 1s shown in Fig. 41. 





Fig 43 —Cambridge Torsiograph (with Cover removed) 


The Cambridge Torsional Vibrograph.—This instrument 
is shown in Figs. 42 and 48. The general working 
principle comprising the light aluminium pulley and the 
floating flywheel, is the same as that of the Geiger, 
but the method of recording is entirely different. 
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The instrument employs, as a recording member, a 
spherical ended stylus tracing on a transparent celluloid 
strip. The principle is based upon the fact that the 
elastic limit of such material as celluloid is low, and 
therefore plastic deformation occurs under comparatively 
shght pressures. The records obtained are so smooth 


; hut ' i ‘7 : W) ‘aaa 





Fig 44 —Cambridge Torsiogram viewed through Mi roscope. 


and fine in character that, when optically magnified, 
readings accurate to about 0-001 millimetre can be made. 

Instruments recording on this principle need make 
only small movements, and consequently the controlling 
forces used in the instrument are relatively great and 
a high frequency response can be obtained, whilst 
inertia reactions are small. 

Fig. 44 shows the appearance of the record when 
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viewed through a simple microscope. The formation of 
the stylus point is important, for if a fine geometrical 
point were used, it would cut or tear a jagged line and 
would require the application of comparatively large 
pressures. It is important that no material be cut or 
torn from the surface. 

For correctly timing the torsiograms, a 1/10th second 
chronometer operates a second stylus, whilst a third 
stylus is operated by contact breaker from the injector 
or exhaust valve. 

To facilitate the examination and enlargement of the 
torsiograms, the portable viewer and enlarger shown 
in Fig. 45 15 used. An enlarged view of the record 
(magnified 10 times) 15 obtained on a focussing sereen 
measuring 5+ < 3} inches. On this screen are ruled 
horizontal lines giving the scale of angular deflection 
in degrees. 

The focussing screen can be replaced by a camcra 
dark slide, cnabling photographic enlargements to be 
obtained on cither a plate negative or on bromide 
paper. From the plate negative, a further enlargement 
can be made if desired, so that records of the original, 
magnified 100 times, can easily be obtained. This is 
particularly valuable when a minute examination of 
the wave form is required. 

A typical torsiogram taken from a diesel engine is 
shown in Fig. 46. 

Operating the Torsiograph.—As in numerous other 
branches of engineering, where cntirely fallacious con- 
clusions may be drawn as a result of incorrect use of 
a measuring instrument, results that are of value can 
be obtained only by knowing, first, the nature of the 
vibration to be measured, and second, the most suitable 
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Fig 45.—Cambridge Portable Viewer and Enlarger 
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Fig. 46.—Typical Cambridge Torsiogram. 
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application of the torsiograph. The general procedure 
can best be set forth as the series of rules given below. 
Some of these rules are incorporated in the instructions 
supplied with the instruments, whilst others are based 
on the author’s expericnee. 

(a) The torsiograph belt should not be provided with 
an ordinary belt fastener, but with a butt joint as 
shown in Fig. £7. 

(b) The material of the belting should be cither the 
stretchless material supplied by the Geiger Co., or steel 
tape 0-003 inch thick, the latter being preferable. 
Ordinary tapes, such as that used for tachometers, should 
not be used. 

(c) Pulley centres should be as short as possible, 


si A a 


Fig. 47.— Butt Jomt Belt Fastener. 


especially in the case of high frequency vibrations. 
The belt should be run sufficiently tightly to obviate 
any possibility of slip. In suitable locations, the 
torsiograph may be direct driven from the end of the 
shaft. 

(d) The ratio of pulley sizes is immaterial. If the 
pulleys differ in diameter, the scale of the torsiogram 
should be multiplied by d,/d,, where d, is the diameter 
of the torsiograph pulley and d, is the diameter of the 
driving pulley. 

(e) The vibrations should, in general, be measured 
at a point of maximum amplitude. This point is 
obtained from the normal elastic curve. 

(f) Where the one and two noded natural frequencies 


are comparatively close together, a diagram taken in 
7 
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compliance with rule (e) will give a curve which is a 
combination of the two modes of vibration and which 
thercfore requires a harmonic analysis to obtain the 
separate records. Under these conditions it is very 
helpful to take additional torsiograms from one of the 
nodes of the two node mode of vibration (see Fig. 48). 
This clears the torsiogram of all but the one node 


vibrations. 





Torsiograph driven trom here does 
not record two node vibrations 


Fig. 48.—Mcthod of Obtaining Clear Record of One Noded Vibrations. 


It oceasionally happens that, when using a torsiograph, 
the stylus vibrates violently between the stops, as 
shown in Fig. 49. The usual cause of this phenomenon 
is the use of too high a magnification, the vibration 
transmitted from the engine sctting up a synchronous 
vibration of the stylus betwcen the stops. 

In cases where a reduction of the magnification does 
not cure the trouble, an examination of the engine 
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will usually reveal the cause. The fault may be found 
in a faulty injector or a sticky or bouncing valve, either 
of which causes, momentarily and periodically reduces 
the power of one cylinder and so causes a sudden speed 
variation. 
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Fig. 49.—Torsiogram showing Violent Instrument Vibration. 


THE INTERPRETATION OF THE TORSIOGRAM 


Havinec obtained a reliable reading from the instrument 
it is now necessary that the torsiogram shall be correctly 
interpreted. In this connection it must be remembered 
that only in critical speeds which are severe compared 
with others in their immediate neighbourhood, is the 
exclusive presence of a particular order of critical 
indicated by the recording of a pure sine wave. In 
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the majority of cases, the diagram of the critical speed 
under investigation has, superimposed on it, the out of 
phase components of neighbouring critical specds. 
This phenomenon produces a distortion in the diagram 
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Fig. 5] —Torstogram from Six Cylinder Engine 


and tends to complicate matters somewhat when 
engine damping factors are being investigated. 

Figs. 50 and 51 are two typical torsiograms. The 
former shows a pure 4th order vibration recorded from 
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a four cylinder engine. The latter, which is recorded 
from a six cylinder engine, shows a 6th order vibration 
m which the out of phase components of the 54th and 
64th order minors produce the distortion mentioned 
above. 

There are two methods of dealing with this latter 
type of torsiogram when an accurate value for the 
engine damping factor is being sought. The first 
consists of analysing the diagram into its constituent 
harmonics, thereby obtaining a series of pure sine 
waves. The second consists of the reverse process. 
By this method a series of torsiograms is taken over 
a given speed range, and a graph of peak amplitudes is 
plotted against the speed. Then, assuming a damping 
factor, the theoretical amplitude curve is built up, 
If this does not agree with the measured curve, a 
further damping factor is assumed until agreement is 
reached. 

The latter method appears at first sight to be the 
more laborious. A number of cases investigated by 
the author, however, show that, if the Shannon method 
of calculation is used, and due allowance is made for 
external damping (if any), exceptionally close agreement 
-between measured and calculated amplitudes is obtained. 
This being so, the harmonic components of any particular 
torsiograi can, if required, be read off from the separate 
components of the theoretical curve and the somewhat 
lengthy calculation called for by harmonic analysis is 
thereby avoided. 

Figs. 52 and 53 show two typical amplitude-speed 
curves, the method of constructing which is as follows :— 

Using the Shannon method of calculation, the 
amplitude at each individual critical speed is calculated. 
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Knowing the dynamic magnifier, Fig. 3 is consulted and 
the amplitude produced by each individual critical over 
the entire speed range is then plotted. These are 
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Fig. 52.—Amplitude-Speed Curve, Six Cylinder Diesel Engine. 
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shown dotted in Figs. 52 and 58. The resultant is then 
found as the algebraic sum of the ordinates. 
In Figs. 52 and 58 the actual measured amplitudes 
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have also been plotted and it will be seen that very 
good agreement has been obtained. 

To obtain the stress-speed curve is now a relatively 
simple matter. From the table for the calculation of 
the natural frequency of the system, the stress S$ per 
degree twist at No. 1 crank is known (see Table I). 
The ordinates of the amplitude-speed curve are thus 
multiplied by S and the required curve is obtained. 

It should be appreciated that in what has been 
said above in connection with the interpretation of 
the torsiogram, the latter has been assumed to have 
been taken from the end of the engine opposite to the 
flywheel. If the torsiogram is taken from elsewhere on 
the system, it is necessary, before applying the above 
analysis, to correct it to the amplitude at No. 1 crank 
by multiplying by the ratio of the ordinates of the 
normal elastic curve at the two points. 


THE ELIMINATION OF TORSIONAL VIBRATION 


THE estimation of probable critical speeds is to-day 
one of the most important aspects of crankshaft design. 
Whilst it may appear a comparatively casy matter first 
to design the crankshaft and then build the remainder 
of the engine round it, circumstances often force the 
hand of the designer and he finds that the crankshaft 
to fit the job has a dangerous major critical on the 
running speed. 

The problem then arises how to modify the crankshaft 
to avoid this critical speed without upsetting the 
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remainder of the design. A number of solutions are 
available, the most important of which are :— 

(a) Alter the crankshaft stiffness. 

(6) Alter the rotating and reciprocating masses. 

(c) Fit a torsional vibration damper. 

In the case of (a) and (6) stiffness should, wherever 
possible, be increased and mass reduced, so that the 
critical speed is raised. If the reverse process is adopted, 
the lowered critical speed has to be run through each 
time the engine is started up and shut down. 

Modifications to crankshaft stiffness usually take the 
form of increased journal diameters. An increase in 
the crankpin diameter, whilst effective, has two dis- 
advantages. It increases the weight of the large end, 
thus partially offsetting the inercase in stiffness, and 
renders more difficult the removal of the connecting 
rod through the cylinder liner. Altcrations to the 
crankwebs are of comparatively little value as increased 
weight again gocs with increased stiffness. 

It should be remembered that alterations to stiffness 
are most effective near the node, the value decreasing 
as the distance from the node increases. 

Modifications to rotating and reciprocating weights 
offer scope only in certain circumstances. Boring the 
crankpins and chamfering the crankwebs are solutions 
which are almost ineffective when there are no balance 
weights, but which are most effective when balance 
weights are fitted. The reason is that in the former 
case the decrease in weight 1s to all intents and purposes 
offset by the reductions in stiffness, whilst in the latter 
case, to preserve the same degree of balance, an 
appreciable amount of matcrial can be removed from 
the outside periphery of the balance weights. 
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Increasing the weight of the flywheel has little effect 
if the moment of inertia is sixty or more times that 
of each cylinder line. In the case of marine engines 
with small flywheels, alterations to the flywheel mass 
are very effective. 

It is not usually desirable to make alterations to the 
reciprocating masses. Apart from the fact that only 
one half of their weight can be assumed to act at 
crankpin radius, considerations such as engine efficiency 
and balance dictate that weight shall be kept at a 
minimum and remedies for critical speeds sought else- 
where. 

Alterations to mass are most effective at points 
remote from the node. Since the node itself is virtually 
an infinite flywheel, the addition or subtraction of mass 
at this point is totally ineffective. 

The Torsional Vibration Damper.—Whilst it is 
axiomatic that alterations to avoid critical speeds 
should be made in the drawing office and not in the 
shops, the designer is often caught unawares when a 
stock engine is sold to run at some speed other than 
standard, or when a normally constant speed engine 
has to run at varying speeds. In either case a serious 
major critical may be encountered, and as alterations 
to mass and stiffness are difficult and sometimes 
impossible, it becomes necessary to fit a torsional 
vibration damper. 

The damper, which was invented by Dr. Lanchester, 
consists essentially of a mass driven through friction 
dises as shown in Fig. 54. The hub a carries the 
driving discs 6, which are free to slide in longitudinal 
grooves in the hub. The driven discs c, which are free 
to slide in longitudinal grooves in the rim d, are pressed 
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against discs b, by means of the springs e. The action of 
the damper is that, should the shaft carrying the hub be 
subjected to violent vibration, slipping occurs between 
the dises b and c, and the vibration is thereby damped 
down by viscous drag. The damper is, of course, always 
ie 
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Fig. 54.—Lanchester Damper. _ Fig. 55.—Solid Friction Damper. 


placed on the shaft at the point of normal maximum 
amplitude. 

The design of a damper from first principles presents 
many difficulties that have yet to be overcome. Major 
Carter has done much research work on the subject 
and has published his results,* and whilst these are 


* B. C. Carter, “ Torsional Vibration in Engines,” Aeronautical Research 
Committee Report No, 1053. 
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admittedly incomplete from the point of view of a 
multi-cylinder diesel engine (the work relates mainly 
to single crank radial engines) certain valuable con- 
clusions have been reached. 

When the friction between the plates is zero, the 
dampcr has no effect and the position and magnitude 
of the critical speeds remain unchanged. As _ the 
frictional resistance increases, the critical speed becomes 
lower and the possible amplitude of the vibrations 
decreases. This effect persists down to a certain point, 
after which a further increase in frictional resistance 
reduces still more the critical speed, but increases the 
possible amplitude of the vibrations, so that when the 
plates are solidly locked together, the damper becomes 
an added flywheel and the possible vibration amplitude 
returns to its original value. 

In order that the damper may be effective, its moment 
of inertia should be at least one half of that of the 
total crank masses. The effectiveness may be increased 
by using a suitable connecting shaft stiffness, or by 
increasing the moment of inertia. 

In theory, the Lanchester damper is of the “ viscous 
friction ”’ type, that is, the friction is independent of 
spring pressure but is proportional to the areca of the 
friction discs. For this reason many of the earlier 
dampers were provided with a large number of discs 
and no means of adjusting the springs. 

Accumulated experience has, however, shown that 
the effectiveness of this type of damper can be varied 
by means of spring pressure. Thus the friction is not 
‘* viscous ”’ but approaches the state known as “‘ solid,”’ 
i.€. independent of the slipping velocity of the rubbing 
members but proportional to the pressure between them. 
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This experience has led designers to search for a 
simpler and cheaper design and the outcome is that 
for medium and high speed engine work, the type 
shown in Fig. 55 1s almost universally adopted. 

In this damper, the friction medium consists of two 
Ferobestos dises 6 which are clamped between the 
hub a and the two dises ec. The dises ¢ are free to rotate 
on the hub, and the latter is firmly keyed to the shaft. 
The pressure of the dises_¢ on the friction medium. 1s 
controlled by the adjustable springs d. 

The action of this damper differs from that of the 
pure viscous friction type in that whereas the latter is 
always slipping so long as there is any vibration, the 
solid friction type acts as an added flywheel until the 
vibration torque exceeds the static torque setting, at 
which point slipping begins. It is thus possible, with 
a constant speed engine having a major critica] at the 
running speed, to fit a solid friction damper sct in such 
a way that it acts as an added flywhcel at normal 
speed, and slips only when starting up and shutting 
down. 

Expressions which greatly facilitate the design of the 
solid friction damper have been evolved by Den Hartog.* 
They are as follows :— 


Energy dissipated per cycle 


1 
éw = git ab? ft. lbs... . (58) 


and the optimum friction torque setting at which this 


dissipation is a maximum 


Tz = Tw ft. Ibs. . =. (59) 


* « Mechanical Vibrations,’ by Den Hartog. 
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where I, = the moment of inertia of the slipping 
portion of the damper in lbs. ft.? sec.? 

a - phase velocity in radians per sec. 
@ — vibration amplitude at the point of applica- 

tion of the damper in radians. 


Assume now that a damper is to be fitted at the 
free end of a diesel engine crankshaft and that the 
stiffness of the connecting shaft is comparatively high, 
i.e. the vibration amplitude at the damper is the same 
as that at No. 1 crank. 

The input work per cyele from the harmonic forces, 
from Equation (46), 1s 

dw =- rH ,,0Xa ft. Ibs. 


The energy dissipated per eyele by the damper is 
1 
= eta ! 2 = e e e 6 
6w, = a0 To?@? ft. Ibs (60) 


If engine damping be neglected, 7.e. if the damper 
is designed to do the whole of the vibration control, 
dw = dw, and 

2 
Li ee Ibs. ft.? sec.? : . (61) 

From the frequency calculation table the stress per 
degree twist at all parts of the shaft being given by the 
last column, a value of 0 is chosen so that the maximum 
stress at the node will not exceed a certain predetermined 
value. By successive approximation Xa, w* and I, can 
now be calculated for any value of H,, t.e. for any order 
of critical speed. 

Once I, and 6 have been fixed, T,, can be calculated 
direct from Equation (59). 

There are two points of particular interest in 
Equation (61) above. The first is that J, is inversely 


THE ELIMINATION OF TORSIONAI, VIBRATION I11 


proportional to 0, which means that to halve the stress 
it is necessary to double the inertia of the damper if 
optimum conditions are to be maintained. The second 
is that J, is proportional to H,2a. This means that 
a dampcr can be designed for one critical speed only. 
For all others, optimum conditions cannot be main- 
tained. 

Although this latter point may appear at first sight 
to present insurmountable difficulties, experienee shows 
that, provided the damper is sct to deal with the most 
severe critical in the running range, no trouble from the 
other criticals need be anticipated. 

The Bibby Detuner.—The Bibby detuner is a device 
which, instead of damping out the torsional vibrations, 
detunes the system in such a way that resonance 
conditions cannot build up. It is, in effect, an ingenious 
adaptation of the well-known Bibby coupling. 

This coupling, which is shown part-sectioned in 
Fig. 56, consists essentially of driving and driven halves 
a and b, each provided with teeth and coupled together 
by means of a continuous spring c, which is wound 
between the teeth in the manner shown. The teeth 
are flared at the ends d in a predetermined manner, 
thus allowing more clearance for the spring at these 
points. 

The principle of operation is that, when no torque 
is being transmitted, the spring takes up its normal 
position with its members parallel to the driving shaft, 
as shown in Fig. 56. Under these conditions, the 
coupling has a torsional stiffness which is dependent 
on the length of spring between the outer ends of the 
teeth e. When overload is being transmitted, the 
driving half of the coupling moves forward axially 
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relative to the driven half, and the spring lies along the 
flared flanks of the teeth for their entire length. In 
this case the torsional stiffness depends on the length 
of spring between the imner ends of the teeth. For 
normal working load, an mtermediate spring position 
and torsional stiffness is used. 





Fig 56 — Bibby Coupling 


The coupling thus possesses the unique property of 
changing its torsional stiffness as the load and consequent 
angular deficction changes. 

The detuncr, in one of its forms, consists of a flywheel 
mounted on the engine crankshaft through a Bibby 
coupling, in the position normally occupied by a damper. 
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The above mentioned property of the spring gives a 
progressive change in the torsional stiffness of the 
system as the amplitude increases. Torsional resonance 
is therefore ‘‘ detuned ” as soon as it attempts to build 
up. The spring characteristic must be designed for the 
particular system to which the detuner is applied, and 
It is possible, by careful design, to detune two and 
sometimes three frequencies. 

In Fig. 57 is shown the amplitude-speed curve for 
the crankshaft of a dicsel engine. The effect of fitting 
a detuner is clearly shown. 
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